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Abstract 

Among the different dynamical processes that take place in polymers, methyl group rotation 
is perhaps the simplest one, since all the relevant interactions on the methyl group can be 
condensed in an effective mean-field one-dimensional potential. Recent experimental neutron 
scattering results have stimulated a new revival of the interest on methyl group dynamics in 
glasses and polymer systems. The existence of quantum rotational tunnelling of methyl groups 
in polymers was expected for a long time but only very recently (1998), these processes have 
been directly observed by high-resolution neutron scattering techniques. This paper revises 
and summarizes the work done on this topic over last ten years by means of neutron scattering. 
It is shown that the results obtained in many chemically and structurally different polymers 
can be consistently described in the whole temperature range — from the quantum tunnelling 
limit to the classical hopping regime — as well as in the librational spectrum, in terms of the 
Rotation Rate Distribution Model (RRDM), which was first proposed in 1994. This model 
introduces a distribution of potential barriers for methyl group rotation, which is associated to 
the disorder present in any structural glass. The molecular and structural origin of the barrier 
distribution in polymers is discussed on the basis of a huge collection of investigations reported 
in the literature, including recent fully atomistic molecular dynamics simulations. 
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I. INTRODUCTION 

The different dynamic processes present in amor- 
phous polymers cover a extremely wide time scale, 
spanning from ca. 10~ 13 s to years (see, e.g., Rcfs. 
[1-4]). These processes include terminal relaxations, 
conformational rearrangements, segmental dynamics, 
localized Johari- Goldstein relaxations, rotations of 
side groups, or fast vibrational dynamics. One of 
the simplest processes among the latters is methyl 
group, -CH 3 , rotation. Many natural and synthetic 
macromolecular chains contain this simple molecular 
unit as a side group or as a part of more compli- 
cated side groups. In the glassy state of polymeric 
materials, at temperatures well below the glass transi- 
tion region, one can assume as a good approximation 
that the main-chain dynamics are completely frozen- 
in and only small side units, as the methyl groups, 
can still remain mobile. In such conditions, the inter- 
action between the methyl group and its environment 
is often well approximated by an effective mean-field 
or "single-particle" rotational potential [5-7]. The 
methyl group can be regarded as a rigid rotor because 
the strength of the C-H covalent bonds allows one to 
neglect the internal degrees of freedom in comparison 
with both translational and rotational motions of the 
group as a whole. Hence, the single-particle potential 
only depends on one characteristic angular coordinate 
4>, which is measured in the plane perpendicular to 
the C3-symmetry axis of the methyl group. The latter 
corresponds to the bond joining the carbon to the rest 
of the molecule (see Fig. 1). 

In the simplest approximation (see Section III) the 
rotational potential is threefold: 



V(0) = ^(l-cos30). 



(1) 



Note that, in this form, V3 corresponds to both the 
maximum and the amplitude of the potential V (</>). 
Figure 1 shows the quantized energy levels of a three- 
fold potential with barrier height V3 = 500 K (in 
the following V3 will be given as V^/fee, with ks the 
Boltzmann constant). Such levels are obtained [5-7] 
by solving the corresponding stationary Schrodinger 
equation, which in the case of a threefold potential 
takes the form of the well-known Mathieu equation 
(see, e.g., Ref. [8]). The energy levels i = 0,1,2... of 
the individual potential wells are named torsional or 
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librational levels, and are split (with energy splitting 
Aj) due to the coupling between the single- well wave- 
functions. Three quantities characterizing the methyl 
group rotational potential are experimentally acces- 
sible (see below): i) the energies E i of the transi- 
tions between the librational levels, specially the tran- 
sition Eqi between the ground and first excited level, 

ii) the energy splitting of the ground librational state, 
Ao = fiujt, with u>t the quantum tunnelling frequency, 

iii) the activation energy Ea for methyl group classical 
hopping between adjacent wells, defined as the difer- 
ence between the top of the barrier and the ground 
state (see Figure 1). The two former quantities can be 
directly accessed by inelastic spectroscopic techniques. 
The latter is indirectly derived from the temperature 
dependence of the hopping rate, which can be deter- 
mined by quasielastic spectroscopy or by relaxation 
techniques (see below). 

According to the scheme shown in Figure 1, methyl 
group dynamics at very low temperature (typically 
T < 20 K) is dominated by quantum rotational tun- 
nelling, whereas classical hopping processes over the 
rotational barrier control the non-vibrational dynam- 
ics at high temperature (typically T > 80 K). Due 
to the one-dimensional character of the rotational po- 
tential, methyl group dynamics is one of the simplest 
examples of barrier penetration phenomena involving 
the interaction of a quantum system with lattice vi- 
brations. For that reason, it has attracted a great 
theoretical interest as an ideal case for the understand- 
ing of the fundamental problem of the transition from 
quantum tunnelling to classical hopping (see, e.g, Refs. 
[9-14]). Moreover, due to the high sensivity of the ro- 
tational potential to the chemical and geometrical na- 
ture of the environment, and to the relative simplic- 
ity of their dynamics, methyl groups can be used as 
built-in probes for exploring structural and dynamic 
properties of the host material. 

Dynamics, and in particular rotational tunnelling, 
of methyl groups has been widely investigated in 
molecular crystals by different techniques over last 
decades, (see, e.g., [5-7]). Measurements of tun- 
nelling frequencies of methyl groups in about 200 sys- 
tems of very different chemical nature have been com- 
piled in Ref. [5]. The used techniques include mostly 
neutron scattering and nuclear magnetic resonance 
(NMR) methods (see, e.g., Ref. [15] for a review on 
NMR measurements), but also calorimetry, electron 



nuclear double resonance, infrared or hole burning op- 
tical spectroscopy. However, no method yields such a 
direct model independent insight as inelastic neutron 
scattering (INS), where the tunnelling frequency is di- 
rectly observed as two resolution-width inelastic peaks 
at ±Twot- Nevertheless, it must be stressed that this 
observation is only possible for moderate and weak po- 
tentials (V3 < 700 K). For stronger potentials Tiuj t is 
smaller than the current available instrumental reso- 
lution (rs 0.3 /tteV) and cannot be detected by neutron 
scattering spectroscopy. In such cases NMR is the 
mostly employed technique. 

The lowest librational transitions E$\, Eq2 take val- 
ues typically in the range 10 — 60 meV and can also be 
directly detected by INS as resolution-width inelastic 
peaks. However, they are observed in the same energy 
range that the other intramolecular vibrational modes 
and the lattice phonons. Though this fact often com- 
plicates the identification of methyl group librational 
transitions, measurements of the latters are useful for 
an accurate determination of the rotational potential, 
since as the tunnelling frequency, they strongly depend 
on the height and shape of the barrier. 

Neutron scattering techniques have also been ap- 
plied over the last decades to investigate methyl group 
dynamics in polymers. It is worth emphasizing that 
polymer systems are in general highly disordered ma- 
terials. This implies that, in principle, one might ex- 
pect the existence of a distribution of potential barri- 
ers for methyl group rotation, resulting from the dif- 
ferent local environments present in the system. In- 
deed, the experimental librational peaks reported in 
the first INS investigations [16,17] in polymers showed 
a broad shape, as might be expected for an inherent 
distribution of potential barriers. However, this as- 
pect of the problem, which had already been consid- 
ered in NMR investigations [18,19], was ignored for a 
long time. Hence, first quasielastic neutron scatter- 
ing data (QENS, see Section II) on classical hopping 
of methyl groups in polymers were analyzed, as for 
the case of crystalline systems, by assuming a unique 
value of the rotational potential. However, this pro- 
cedure provided, in contrast to the result observed in 
crystalline systems, a non-Arrhenius temperature de- 
pendence of the hopping rate, and an apparent tem- 
perature dependence of the geometry of the motion, 
with large deviations from threefold rotation [20-26]. 

The presence of quantum tunnelling effects at un- 
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usually high temperature was considered as an unreli- 
able explanation of the observed features, and a pos- 
sible interpretation, by mantaining the picture of a 
unique rotational barrier, was discussed in terms of 
sixfold symmetry contributions to the rotational po- 
tential [20,21]. Another explanation was introduced 
in terms of temperature dependent fractions of rotat- 
ing and non-rotating methyl groups [22,27]. Though 
a reasonable description of the experimental data was 
achieved, a physical justification of the latter hypoth- 
esis was lacking. Finally, it became clear that the 
basic assumption of a unique barrier was wrong for 
amorphous polymers when experiments on a same sys- 
tem performed in different spectrometers, i.e., in dif- 
ferent time/energy windows, provided uncompatible 
results for the measured hopping rate [23,24]. In the 
seminal works of Refs. [25] and [26], QENS data for 
methyl group hopping in polymers were successfully 
interpreted by introducing a broad distribution of ro- 
tational barriers, resulting from the highly disordered 
nature of the glassy state. In last years, a noticeable 
amount of experiments in different systems have sup- 
ported this picture (see Sections V and VI). 

Another aspect of the problem of methyl group dy- 
namics in polymers is the dynamic behavior at very 
low temperature. From the early times, there was in 
the literature a number of claims or predictions about 
low temperature quantum effects. In particular, NMR 
measurements by Hoch et al. [28] showed methyl group 
dynamics in poly (vinyl acetate), PVAc, to be active 
down to very low temperature. Rotational tunnelling 
of methyl groups was also claimed to be present in 
poly (methyl methacrylate), PMMA, which would ac- 
count for the reported non-Arrhenius behavior of the 
mechanical relaxation time at very low temperature 
[29,30]. However, a quantitative and unambiguous 
interpretation of the experimental data was not pro- 
vided. As mentioned above, the rotational tunnelling 
picture has been evidenced in a vaste collection of crys- 
talline systems by different techniques, and specially 
by neutron scattering, where the tunnelling frequency 
is directly observed as two resolution-width inelastic 
peaks at ±hcj t . For the case of polymers, one might 
expect to observe broad peaks as a consequence of the 
distribution effects. However, observation of peaks, or 
of any scattered intensity over the instrumental reso- 
lution at very low temperature was not reported for a 
long time. 



Finally, high resolution measurements in PVAc [31] 
and PMMA [32] at T w 1 K showed, instead of well- 
defined peaks, a broad scattered intensity similar to 
that commonly observed in the high temperature hop- 
ping regime. This result, which might have been inter- 
preted as a signature of hopping events at unusually 
low temperature, was indeed explained in terms of ro- 
tational tunnelling, namely as the result of the pres- 
ence of a strongly assymetric distribution of tunnelling 
frequencies, with the maximum located well beneath 
the instrumental resolution. The latter distribution 
just followed from the distribution of potential barri- 
ers for methyl group rotation. This interpretation was 
supported by exploting the expected isotope effect on 
the tunnelling frequencies [32] (see Section VI). The 
consistency of the distribution picture, formalized in 
the Rotation Rate Distribution Model (RRDM, see 
Section IV) was supported by the fact that the same 
barrier distribution was able to quantitatively account 
for the experimental spectra in all the temperature 
range [33,34], covering the crossover from rotational 
tunnelling to classical hopping, as well as the libra- 
tional peaks observed by INS. 

The grounds of the RRDM can, in principle, be ex- 
tended to non-polymeric disordered systems, as has 
been successfully checked in the glassy state of toluene 
[35] and sodium acetate trihydrate [36] (see Section 
VII). Contrary to the case of polymer systems, results 
for methyl group dynamics in molecular glasses can be 
directly compared with those for the crystalline state, 
allowing one to get insight into the molecular origin of 
the barrier distribution. 

The goal of this review is to summarize the progress 
in the topic of methyl group dynamics in polymers 
that has been achieved over last 10 years by means of 
neutron scattering techniques, together with the de- 
velopment of the RRDM. The molecular origin of the 
distribution of potential barriers, which is an essential 
ingredient of the RRDM, is discussed on the basis of 
the experimental and computational investigations re- 
ported in the literature. A brief introduction to the 
theoretical and experimental grounds of neutron scat- 
tering techniques is also given in next Section. 
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II. NEUTRON SCATTERING TECHNIQUES 
AND INSTRUMENTATION 



A. Neutron scattering techniques 

Neutrons are adequate probes for the investigation 
of both structure and dynamics of condensed matter. 
There are three major features for that: i) their typ- 
ical wavelengths are of a few A, which correspond to 
the typical interatomic distance, ii) their typical ener- 
gies are of the order of some meV, which correspond 
to the energy scale of the typical excitations in con- 
densed matter, iii) neutrons being particles without 
electrical charge, they can penetrate the sample and 
provide information about bulk properties, in contrast 
to charged particles like electrons, which mainly probe 
surface properties. 

The range of the nucleus-neutron interaction is ~ 
1.5 x 10 ~ 5 A, much smaller than the neutron wave- 
length and the nuclear radius. Hence, scattering can 
be approximated as isotropic, and the interaction is 
modelled by the Fermi pseudopotential: 



V F 



2nfr 



W(r-R), 



(2) 



where m and r are respectively the neutron mass and 
position, and R is the nucleus position. The scattering 
length operator, 6, is given by 



6=6. 



2h, 



coh 



si, 



(3) 



where S and I are respectively the neutron and nuclear 
spins, and b co h and bi nc are respectively the coherent 
and incoherent scattering lengths. 

As in any scattering technique, there are two ba- 
sic quantities to be measured by neutron scattering: 
the energy transfer, hu> — ti 2 (k 2 — fc§)/2m, which is 
the difference between the final and the initial neutron 
energy, and the momentum transfer, ?iQ = Hk — Tik , 
where k and k are respectively the final and initial 
neutron wavevectors (see Fig. 2). The number of neu- 
trons scattered within a solid angle between ft and 
+ dSl, which have changed their energy in hcu and 
their momentum in HQ, is proportional to the double- 
differential cross-section [37-40]. The latter can be 
split in two parts, the coherent and the incoherent 
double-differential cross-sections: 



which can be expressed as a function of the scattering 
functions: 



8 2 a \ 
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(5) 



In this equation a co h and cr inc are respectively the co- 
herent and incoherent cross-sections. S co h{Q,w) and 
Sinc(Q,u) are respectively the coherent and incoher- 
ent scattering functions. These functions give account 
for the amplitude transitions, mediated by the Fermi 
pseudopotential, between the eigenstates of the sam- 
ple. While the coherent contribution results from in- 
terferences between the neutron waves dispersed by 
the nuclei, the incoherent contribution is equivalent 
to the dispersion of the isolated nuclei. In the clas- 
sical limit, S co h{Q,u) and Si nc (Q,oj) are related via 
double Fourier transform, respectively, with the pair- 
correlation and self-correlation functions of the atomic 
positions. 

Protium hydrogen (H) has a large incoherent scat- 
tering cross-section as compared to its coherent part 
and to the total cross-section, a tot — ai nc + cr co h, of the 
other elements typically present in polymer materials 
(see Table I). Hence, a neutron scattering experiment 
on a protonated polymer sample provides mainly in- 
formation on the self-correlation function of the pro- 
tium hydrogens. Moreover, due to the large difference 
between the incoherent cross-section of protium and 
the total cross-section of deuterium (D), the intensity 
scattered by a subset of the protium hydrogens can be 
strongly enhanced if selective deuteration of the other 
hydrogens is possible by means of chemical methods. 

When neutron scattering is used for dynamic inves- 
tigations, one finds three different types of scattering 
processes: elastic, quasielastic and inelastic scatter- 
ing. Inelastic neutron scattering (INS) corresponds 
to the case where the neutron exchanges energy with 
the sample excitations, and manifests in the experi- 
mental spectrum as two resolution-width lines of fi- 
nite energy (for neutron energy gain and loss). On 
the contrary, elastic scattering corresponds to the case 
where the neutron energy remains constant, yielding a 
resolution-width elastic line. Finally, quasielastic neu- 
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tron scattering (QENS) is that associated with the en- 
ergy transfer resulting from the Doppler effect occur- 
ring when a neutron interacts elastically with a moving 
particle, and manifests as a broad line around the elas- 
tic line. Of course, inelastic or quasielastic processes 
are no detectable when the involved energy change 
is smaller than the instrumental energy resolution, 
yielding an additional contribution to the experimen- 
tal elastic line. Quasielastic and inelastic processes 
typically cover a range \Uw\ < 2 mcV and \Twj\ > 2 
meV respectively, the latter corresponding to typi- 
cal vibrational excitations in condensed matter. The 
terms QENS and INS are often used to refer to these 
energy ranges and the dynamic window of the corre- 
sponding spectrometers. However, it must be stressed 
that excitacions in the ^eV range can also be present 
in the system, as the case of methyl group tunnelling. 
Hence, such excitations, which strictly correspond to 
INS, are accessible by means of instruments specifi- 
cally designed for QENS. 

In the liquid state, molecules are able to diffuse and 
loss memory from their initial positions. Molecular 
diffusion yields no elastic component in neutron scat- 
tering spectra. On the contrary, if the motion of a 
group of particles is constrained into a localized region, 
as the case of methyl group rotation in a crystalline or 
glassy system, in addition to the eventual inelastic and 
quasielastic contributions, such a motion provides an 
elastic componet in the experimental spectrum. The 
clastic peak is modulated by a Q-dependent quantity 
named the elastic- incoherent-structure- factor (EISF), 
A(Q). The EISF gives account, via Fourier transform, 
for the average positions of the individual particles 
performing the localized motion, i.e, it provides infor- 
mation about the geometry of the motion. 

Integration of the coherent scattering function pro- 
vides the elastic coherent structure factor, S(Q). 
Fourier transformation of the latter gives account for 
the probability distribution of the configurations of all 
the ensemble of particles in the system, i.e., it provides 
static structural information. 

B. Neutron scattering instrumentation 

Neutron scattering investigations on methyl group 
dynamics in polymers are carried out mainly by means 
of time-of-flight (TOF) and high-resolution backscat- 



tcring (BS) techniques. In addition, in most of the 
cases it is convenient to evaluate the total coherent 
and incoherent contributions to the scattered intensity. 
For this purpose, diffraction experiments with spin po- 
larization analysis are neccessary. The showcase pre- 
sented here for PVAc (see Section V) corresponds to 
experiments performed at the spectrometers IN5, IN6, 
IN16 and D7 of the Institute Laue Langevin, (ILL, 
Grenoble, France) and TOSCA at the source ISIS of 
the Rutherford Applcton Laboratory (RAL, Chilton, 
United Kingdom). Next we summarize the technical 
characteristics of these spectrometers. A detailed de- 
scription can be found in Refs. [41-46]. 

IN5 [41] and IN6 [42] are TOF spectrometers. By 
means of a monochromator, different wavelengths, A, 
can be selected in the ranges 1.8-16 A for IN5and 4-6 
A for IN6. Most of the data presented here were ob- 
tained with A = 6.5 A at IN5 and A = 5.1 A at IN6. In 
this situation the elastic energy resolution (half-width 
at half maximum, HWHM) is w 25 /xeV for IN5 and 
i=a 50 fieV for IN6. The explored energy window (in 
neutron energy gain) is -0.6 meV < Tuo < 19.5 meV for 
IN5, and -2 meV < hu < 1600 meV for IN6. The scat- 
tering function S(Q,u>) at constant Q can be obtained 
from the experimental quantity which is directly mea- 
sured, namely S , ($, u>), where $ is the scattering angle, 
by an interpolation procedure according to the trivial 
geometrical relation: 

fi 2 Q 2 /2m = 2E + huj- 2^/E (E + ?iw)cos$, (6) 

where m is the neutron mass and E = h 2 kl/2 rn 
is the incoming neutron energy The interpolation 
procedure provides reliable spectra in the quasielas- 
tic energy range \huj\ < 2 mcV for a Q-range of 
0.6A < Q < 1.8 A , both ranges being of inter- 
est for investigations on methyl group rotation. 

IN5 and IN6 arc suitable to investigate methyl group 
rotational hopping at intermediate and high temper- 
ature, where the corresponding time scale is typically 
of the order of 10 ps, but usually are not adequate for 
the investigation of rotational tunnelling. Only when 
the rotational barrier is very low (V3 < 200K) the tun- 
nelling lines are accessible in the energy range covered 
by TOF spectrometers. For higher barriers, the use of 
backscattering instruments is needed. 

IN16 is a high-resolution BS spectrometer [43] which 
can explore excitations in the ^.eV-range, and corre- 
spondingly, time scales of the order of 1 ns. The high 
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resolution of IN 16 makes it adequate for the investiga- 
tion of methyl group dynamics in, both, the classical 
hopping regime at intermediate temperature, and the 
quantum tunnelling regime at very low temperature. 
By using a Si (111) monochromator the wavelength 
of the incident neutrons is A = 6.27 A. The corre- 
sponding accessible Q-range is 0.2 - 1.9 A . When 
deformed Si (111) crystals are used as analyzers a res- 
olution of nearly Gaussian shape with HWHM w 0.45 
/ieV is attained. In the standard operation mode the 
monochromator is subject to an oscillatory motion at 
a frequency of s=s 14 Hz, which allows one to explore an 
energy transfer range up to 15 ^eV. As in backscatter- 
ing spectrometers the incoming beam energy is much 
larger than the maximum energy transfer, Eq. (6) can 
be approximated by a direct relation between the posi- 
tion of the detector and the momentum transfer, given 
by: 



Q = (4tt/A) sin($/2). 



(7) 



BS spectrometers as IN16 can also be used in a differ- 
ent operation mode, known as the fixed-elastic-window 
(FEW) technique. In this technique the monochroma- 
tor is held fixed and the sample temperature is varied 
at a moderate rate. In this way one measures the elas- 
tic intensity as a function of the temperature, since 
the FEW technique only detects the fraction of neu- 
trons with an energy exchange within the instrumental 
resolution. Hence, what is considered as elastic inten- 
sity depends on the resolution of the used spectrom- 
eter. For IN16, with a resolution function of HWHM 
w 0.45 ^eV, the dynamics of scattering particles be- 
ing slower than a characteristic time r w 10 ns are 
observed as elastic scattering. When increasing tem- 
perature, T, the fraction of particles slower than r 
decreases, yielding a decrease of the elastic intensity. 
This decrease is exponential if dynamics in that time 
scale is purely vibrational, and it is given by the in- 
dependent Debye- Waller factor e~ 2W( ^ [37-40], with 
W(Q) = Q 2 (u 2 )/3. In this expression (u 2 ) is the 
mean-squared displacement per particle. For vibra- 
tional motions (u 2 ) oc T. The presence of other dy- 
namic processes is evidenced by deviations from ex- 
ponential behavior, yielding a step-like decay of the 
elastic intensity (see Section VI). 

The spectrometer D7 combines TOF methods with 
spin polarization analysis capabilities [44,45]. Among 
the different types of experiments which are possible 



with D7, that of particular interest for the investiga- 
tion of methyl group dynamics in polymers is the mea- 
surement of the coherent and incoherent contributions 
to the total scattered intensity (without energy anal- 
ysis). These measurements are usually performed at 
low temperature to ensure that no significant intensity 
is lost due to the finite range of energy transfer cov- 
ered by the instrument. By using a neutron incident 
wavelength of 4.8 A, the relevant accessible Q-range, 
Q < 2.5A , corresponds to that covered in dynamic 
experiments with TOF and BS instruments. 

In TOF instruments as IN5 or IN6, high values of 
vibrational frequencies are only accessible in neutron 
energy gain. Since the thermal statistical weight of 
high-energy vibrational excitations is negligible at low 
temperature, such excitations can only be detected by 
the former instruments at high temperature. This be- 
comes a problem for the observation of methyl group 
librational peaks, which, due to anharmonicity effects, 
show a strongly damped intensity at high temperature 
[47-50]. This problem can be avoided by using an in- 
strument which has access to vibrational frequencies 
in neutron energy loss, i.e., which is able to create 
the corresponding vibrational excitations in the sam- 
ple even at very low temperature. An example of a 
TOF spectrometer that is optimized for vibrational 
spectroscopy is TOSCA [46], which allows one to mea- 
sure neutron energy loss in a wide range of values, typ- 
ically 1 meV < huj < 250 meV. TOSCA design is such 
that there is a unique value of Q for each value of the 
energy transfer, so all the detectors are summed up to 
obtain the vibrational spectrum. The arrangement of 
the spectrometer is such that the neutrons are detected 
at a fixed scattering angle $ = 135°. In this situation 
Eq. (6) reduces approximately to the relation Q 2 oc ui 
for low values of the neutron energy transfer (hu> < 40 
meV). Hence, TOSCA arrangement provides a nearly 
direct access to the generalized vibrational density of 
states (VDOS), Z(w), given by [37-40]: 



S(Q,lo)uj 
Q 2 [n(w) + 1]' 



(8) 



where n(u>) = \exp{TvjO /ksT) — 1] _1 is the statistical 
Bose occupation factor. It must stressed that the gen- 
eralized VDOS Z(uj) is not proportional to the real 
VDOS, but a "deformed" function, where excitations 
are weighted by scattering cross-sections, molecular 
masses, and coordinates of the normal modes. Only 



8 



for a strictly monoatomic sample Z(ui) is proportional 
to the real VDOS. In the general case, though the 
experimental peaks will correspond to the real excita- 
tions, their amplitudes will be controlled by the men- 
tioned factors, being larger for the modes involving 
motions of protium hydrogens. 

C. Final experimental considerations 

Instrumental resolution functions are usually de- 
termined from measurements on a vanadium sam- 
ple, which provides solely elastic scattering. The 
same measurements can be used to calibrate detec- 
tors. In order to obtain the actual experimental scat- 
tering function S(Q,ui), raw data must be corrected 
for detector efficiency, sample container and absorp- 
tion. Due to its small scattering cross-section, alu- 
minium is the usual material for sample holders. Slab 
geometry is usually preferred. The typical area of the 
neutron beam is of 10— 15 cm 2 . The sample thickness 
is selected to provide a transmission of 90 — 95 %. For 
fully protonated polymers this condition is attained 
with sample thickness of ~ 0.1 mm. High transmis- 
sion allows one to neglect multiple scattering effects, 
provided that the very low Q-range is not considered 
in the data analysis. This is the case for methyl group 
dynamics, since the EISF at low Q is very close to 
unity (see Fig. 3), and correspondingly, the quasielas- 
tic intensity — which is weighted by 1 — A(Q), see 
Section III — is very weak and difficult to analyze. 

III. METHYL GROUP DYNAMICS IN 
CRYSTALLINE SYSTEMS: THEORY AND 
SCATTERING FUNCTIONS 

A. Rotational tunnelling 

As mentioned in the Introduction, methyl group mo- 
tion at very low temperature (T ~ 1 K) can be mod- 
elled as that of a rigid rotor in an effective single- 
particle one-dimensional rotational potential V(4>). 
The corresponding Hamiltonian is given by [5-7] : 

Hk = -B^ + V{4>). (9) 

In this equation B = ft 2 /2I is the rotational constant 
of the rigid rotor, with / the moment of inertia of the 



methyl group around its threefold symmetry axis. For 
a fully protonated methyl group B = 0.655 meV. For 
a fully deuterated one, B = 0.328 meV. 

The rotational potential V(cj>) must be invariant un- 
der 2n7r/3 rotations, with n — ±1,±2, ±3..., because 
its symmetry cannot be lower than that of the rigid 
rotor, and moreover, both symmetries cannot be un- 
compatible. Hence, V{<f>) can be expanded as a Fourier 
series of 3n-fold terms of amplitude Vs n : 

oo 

v^)^j2^ l - cos ( 3n ^+ s ^ ( 10 ) 

n =l 1 

where 5^ n are constant angular offsets. Indeed, it can 
be formally demonstrated that the rotor-lattice inter- 
action can be expressed in this way provided that V{4>) 
is given by a sum of two-body additive potentials [51]. 
As mentioned in the Introduction, in a first approx- 
imation only the threefold term of the expansion is 
retained: 

If higher-order terms are neccesary to account for 
the experimental results, they often contribute just as 
small corrections to the main threefold term [5]. 

As mentioned in the Introduction, the coupling of 
the wavefunctions of the individual wells splits the li- 
brational levels. The resulting sub-levels are denoted 
as A and E, the latter consisting of the degenerate 
doublet (E a , E^). These labels correspond to the ir- 
reducible representations of the symmetry group C3 
of the Hamiltonian Hr. The eigenstate A remains 
unchanged under a rotation 27r/3, while such an op- 
eration adds a phase ±27r/3 to the spatial part of the 
wave functions of E a and Eb respectively. An eigen- 
state of Hr is characterized by two discrete quantum 
numbers, namely the librational, n, and symmetry, S, 
indexes: \nS). As shown in Fig. 1, the levels A and E 
alternate between consecutive librational levels. 

At very low temperature, T~ IK, fc^T <C E ni and 
the occupation probability of the excited librational 
states i — 1,2,3... is negligible. Hence, the rotor os- 
cillates coherently between the three wells with an os- 
cillation frequency equal to the tunnelling frequency 
uj t . If the latter takes a value between the highest res- 
olution available in BS spectrometers, hui t ~ 0.3 /xeV, 
and the free rotor limit, huj t — B, it is detected as two 
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resolution-width inelastic peaks at oj = ±u> t (see Fig. 
4a). 

The incoherent scattering functions for rotational 
tunnelling in different symmetries can be calculated 
by using symmetry adapted scattering operators [6]. 
Hence, the incoherent scattering operator 



W 



N 

£ 



(S-I^expHQ-r.) 



(12) 



where rj, Ij are the hydrogen positions and spins, and 
S is the neutron spin, can be rewritten, for threefold 
symmetry, as: 



W = W A + W Ea +W Eb . 



(13) 



The eigenstates \A), \E a ) and \E b ) are unchanged by 
Wa , and are permutated cyclically by W Ea , and anti- 
cyclically by W Eb . Hence, the decomposition (13) al- 
lows for a straightforward calculation of the VT-matrix 
elements between the different eigenstates [6]. After 
averaging over all the orientations of the scattering 
vector Q, the following transition amplitudes are ob- 
tained: 

i) f AE {Q) = [l-io(Qr)]/3 for A <- E {E a ,E b ) transi- 
tions. The corresponding contribution to the scatter- 
ing function is 2fA. E (Q)5(u> ± u>t), since there are two 
inelastic transitions lo = +uj t and other two u> = —u>t- 

ii) f EaEb (Q) = [1 - MQr)}/3 for E a <-» E b tran- 
sitions. These two elastic transitions contribute as 
2f EaEb (Q)5(iu). 

hi) fNSc(Q) = [1 + 2jo(<3r)]/3 for transitions without 
symmetry change. These three elastic transitions con- 
tribute as ^fNSc(Q)S( UJ )- 

In these expressions, jo{Qr) — sin(Qr)/Qr, with r the 
H-H distance in the methyl group, (r = 1.78 A). Af- 
ter normalization to one-hydrogen scattering, the cor- 
responding incoherent scattering function for methyl 
group tunnelling is [5-7]: 



(Q,w) 



5 + 4j (Qr) 



S(uj) + 



2[1- jo (Qr) 
9 



[5(w + ut) +S(oj -u> t )}. 



(14) 



The EISF, A(Q) = [5+4j (Qr)]/9, and [l-A{Q)}/2 = 
2[1 — jo{Q r )]/9, which modulate respectively the elas- 
tic and inelastic contributions, are shown in Fig. 3. 



B. Crossover from rotational tunnelling to 
classical hopping 

The usual temperature dependence of neutron scat- 
tering spectra for methyl group dynamics in crystalline 
systems is illustrated in Fig. 4 for crystalline sodium 
acetate trihydrate [36]. At T ~ 1 K (Fig. 4a) the 
spectrum shows a central elastic peak and two inelas- 
tic peaks of resolution width, i.e., delta functions con- 
voluted with the instrumental resolution. According 
to the picture exposed above, the inelastic peaks corre- 
spond to the two transitions A <-» E within the ground 
librational state. The elastic peak includes the ground 
state transitions E a <-» E b and the transitions without 
symmetry change. It also includes the elastic inco- 
herent contribution from other atoms different from 
methyl group hydrogens, and the total coherent con- 
tribution (see below). 

When increasing temperature [36,48,50,52-59], the 
inelastic peaks progressively broaden and shift to the 
clastic peak. At the same time, a quasielastic com- 
ponent appears around the elastic peak, and broad- 
ens with increasing temperature (Figs. 4b, c). Finally, 
in a narrow temperature interval (AT ~ 7 K), the 
quasielastic and inelastic components merge in a sin- 
gle quasielastic line (Fig. 4d). In the following the 
width of this interval will be neglected and repre- 
sented as a unique temperature T c , that will be de- 
noted as the crossover temperature. Above T c , a clas- 
sical picture of thermally activated hopping over the 
rotational barrier is commonly accpted [6,37]. The 
observed single quasielastic line is well described by a 
Lorentzian function. The corresponding HWHM fol- 
lows an Arrhcnius-like temperature dependence: 



r = r oc cM-E A /k B T). 



(15) 



where the classical activation energy Ea , as mentioned 
in the Introduction, is defined as the difference be- 
tween the top of the barrier and the ground state (see 
Fig. 1). In the following Ea will be given in units 
of k E - The preexponcntial factor Too is temperature 
independent and typically takes values of ~ 5 — 10 
meV. 

While rotational tunnelling and classical hopping 
are well established pictures respectively in the low 
and high temperature limits, no satisfactory theory 
for the crossover regime between both dymanic limits 
has won general acceptance [9-14,60-65] . Every model 
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for the crossover regime must give account for the fol- 
lowing phenomenology observed in a large amount of 
experiments [36,48,50,52-59]: 

i) Broadening of the quasielastic and inelastic 
crossover components are well described by Lorentzian 
functions of the same intensity. The relative in- 
tensities to the elastic peak are the same than 
in the purely rotational tunnelling regime, i.e., 
2/a_e(<5)/3/atsc( ( 3) for each of the inelastic lines and 
2fE a E b {Q)/3fNSc{Q) for the quasielastic line. The 
temperature dependence of the Lorentzian HWHM 
follows Arrhenius-like behavior. The observed acti- 
vation energy corresponds approximately to the first 
librational energy: 

r, = 7j exp(-£ i/fc B T). (16) 

T q = lq cxp(-E 01 /k B T). (17) 

In these equations Ti and T q are respectively the 
Lorentzian HWHM for the inelastic and quasielastic 
components. The preexponential factors ji and j q are 
temperature independent, and typically take values in 
the range 0.05-1 meV, much lower than the observed 
values for the classical (see above). It is usually 
found that 7 9 < 7», though both factors take values of 
the same order of magnitude (typically "f q /ji ~ 0.5). 

ii) Shift of the tunnelling peaks also follows Arrhenius- 
like behavior: 

hAut = is exp(-E s /k B T), (18) 

where Au) t is the shift (in absolute value) of the tun- 
nelling peak and j s a temperature independent preex- 
ponential factor, which takes values similar to 7$ and 
7q. The activation energy Es is generally smaller than 
E i though it takes a close value, typically Es/E i ~ 
0.7. Fig. 5 shows all these features for crystalline 
sodium acetate trihydrate as a standard showcase [36] . 

All the theoretical models that aim to give account 
for these universal features state that its physical ori- 
gin is the coupling of the methyl group to the lattice 
phonons, which progressively breaks the coherence of 
the wavefunction of the rotor, leading to incoherent 
tunnelling [9-14,51,60-84]. The standard approach, 



and the only one that is able to give account quantita- 
tively for the observed behavior, at least below T c , is 
to describe the system rotor-lattice by a Hamiltonian: 

H = H R + H B +H C . (19) 

In this equation Hr is the rigid rotor Hamilto- 
nian (9), H B is a bath of harmonic oscillators 
[9,11-13,64,71-73,75-79,82,84]: 

N ( h 2 d 2 1 \ 
ff B = E(-2^^ + 2 m ^J' (20) 

and He is a coupling term between the rotor and the 
bath: 

N 

H c = ^gkXk cos(3</> + 5 k ). (21) 
fe=i 

In these equations Xk, nik, w fc, 9k and Sk are respec- 
tively the displacements, masses and frequencies of the 
harmonic oscillators, the coupling constants and the 
angular offsets relative to the rotor. 

The total Hamiltonian (19) retains the threefold 
symmetry of the interactions methyl group - envi- 
ronment. Hence, an argumentation on the basis of 
transitions between different symmetry states is still 
possible. The eigenvectors of the Hamiltonian (19) 
are characterized by a large set of quantum numbers. 
However, since H is, as Hr, unvariant under the sym- 
metry group C3, the symmetry indexes A, E a and Eb 
are still good quantum numbers for H, and the sym- 
metry adapted formalism can still be used to calculate 
the scattering function. In the limit N — > 00, H has a 
continuous number of eigenstates. Hence, transitions 
with symmetry change will take place not between two 
discrete eigenstates, but between two continuous sets 
of eigenstates. The delta functions for the low temper- 
ature limit [Eq. (14)], corresponding to transitions be- 
tween two eigenstates characterized by a symmetry in- 
dex, are now substituted by Lorentzians, correspond- 
ing to transitions between two large sets of eigenstates, 
each set characterized by a symmetry index, i.e., 

m -+ \es'), (22) 

where S and S' are the initial and final symmetry 
states, and £ and £' a large set of quantum numbers 
characterizing the initial and final state of the bath. 



11 



Within this picture, the inelastic Lorentzians involve 
transitions A <-» E, while the quasileastic Lorentzian 
involve transitions E a <-> Hence, Eqs. (16,17) can 
be rewritten as: 



^ae = jAECVp(-E 01 /k B T) 



?E a E b = lE a E b CXp(~E 01 /k B T). 



(23) 



(24) 



It can be demonstrated that elastic transitions with- 
out symmetry change do not lead to broadening [13], 
yielding an elastic line as in the low temperature limit. 
Taking all this arguments into account, the incoherent 
scattering function for the crossover regime is given by 
[48,77]: 

[L(u; r EaEb ) + L(u + (wt- Auj t )-T AE ) + 

L(co-(ujt-Awt);TAE)}, (25) 

where L(io—v\ rj) is a Lorentzian function of area unity, 
centered at lo = v and HWHM = 77, i.e.: 



1 



L(uj - 1/577) = - 



V 



7r (u) — v) 2 + rj 2 



(26) 



It must be stressed that the functional form (25) is 
a direct consequence of the threefold symmetry of 
the Hamiltonian (19). It is the description of the 
temperature dependence of shift and broadening of 
the Lorentzian lines which depends on the particu- 
lar model for the coupling of the rotor to the bath 
[Eqs. (20, 21)]. In the different models based on this 
approach, Eqs. (23, 24) are obtained as effective ex- 
pressions from a sum over complicated functions con- 
taining Bose factors resonant on the librational fre- 
quencies, the density of states Z(u>) and the coupling 
constants {gk}' 



N 



r ee =^2^'[Z(E j),{gk},n(E i 



0j) 



(27) 



with ££' = AE or E a Eb, and n(ui) the Bose occu- 
pation factor. In practice, for the temperature range 
T < T c , where the crossover regime is experimentally 
observed, one has fc^T <C Eq\. Hence, the sum (27) 
is dominated by the term in E01 and can be approx- 
imated as T^> = 7^/ exp(— Eoi/ksT), in agreement 
with Eqs (23, 24). 



For the case of the shift of the tunnelling frequency, 
the sum is done over all the phonon spectrum: 

/•OO 

hA(j t = dujg s [Z(Lo),{g k },n(uj)}. (28) 
Jo 

Again for k B T <C £01 the integral can be approx- 
imated as an effective Arrhenius equation as (18). 
Since the integration is done over all the phonon spec- 
trum, it is dominated by low-energy contributions, 
yielding an effective activation energy E$ smaller but 
close to £01. 

There is a strong controversy about the validity of 
the model given by Eqs. (19,20,21) for arbitrarcly high 
temperatures. Some authors have claimed [64,77] that 
it is a valid picture even above T c , the temperature 
above which spectra are well described by the classi- 
cal hopping model. According to this statement, the 
observed classical-like Arrhenius behavior (15) is also 
obtained as an effective expression of the sum (27), 
where the highest order contributions would lead, at 
high temperature, to the observed apparent activation 
energy Ea, which is higher but close to the highest li- 
brational transition (see Fig 1). 



C. Classical hopping 

The incoherent scattering function above T c can be 
derived in a classical hopping picture, disconnected 
from the low temperature quantum features above de- 
scribed. A master equation [6,37] is assumed for the 
classical probabilities (denoted as p\, P2, P3) of find- 
ing the rotor in each of the three equilibrium positions 
0,±2tt/3: 



Pi = 


1 

T 


-Pi + 


P2 

2 


+ ?] 


P2 = 


1 

T 


-P2 + 


P3 

2 


+ t! 


P3 = 


1 

r 


~P3 + 


Pi 
2 


+ f] 



(29) 



The hopping time between two equilibrium positions 
is assumed to be negligible in comparison with the 
residence time, t, between consecutive jumps. This 
assumption is confirmed (see e.g., [85-88]) by molec- 
ular dynamics simulations (MDS), which allow one to 
observe the motion of each individual methyl group in 
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the system. Fig. 6 shows computed angular positions 
for a methyl group hindered by a moderate barrier. As 
can bee seen, sudden jumps occur between the three 
equilibrium positions, with long intervals between con- 
secutive jumps. It is worth emphasizing that r is an 
average value over all the time intervals between con- 
secutive jumps of a same methyl group. The fact that 
the latters take different values is a consequence of the 
stochastic nature of the hopping process. For a crys- 
talline system with a unique barrier, r will be the same 
for all the individual methyl groups. For a disordered 
system, r will be different for each individual methyl 
group, as will be discussed in Section IV. 

After solving the master equation (29), averaging 
over all the orientations of Q, and Fourier transform- 
ing to the frequency domain, the following expression 
is obtained [6,37]: 



2[l-j Q (Qr)} 



L(u;;T), 



(30) 



where T — 3t/2 is just assumed to follow the classi- 
cal Arrhenius-like behavior of Eq. (15). The EISF, 
A(Q) = [1 + 2j (Qr)]/3, and 1 - A(Q) = 2[1 - 
3o(Q r )]/3, which modulate respectively the elastic and 
quasielastic contributions, are shown in Fig. 3. 



D. Additional remarks 

It must be noted that, as required by a consis- 
tent description of the temperature dependence of the 
spectra, the scattering function (25) for the crossover 
regime is reduced to Eqs. (14) and (30) respectively 
in the low and high temperature limits. For very low 
temperature, shift and broadening in Eqs. (18,23,24) 
tend to zero, recovering Eq. (14). On the other hand, 
tunnelling is supressed above the crossover tempera- 
ture, i.e., Aw t = uj t , and the three Lorentzians merge 
into a single line if T ae = YE a E b = T. Hence, Eq.(30) 
is recovered. 

Finally, the total scattering function S{Q,u) is ob- 
tained by adding the incoherent contribution of the 
other atoms, and the total coherent contribution. If, 
except for methyl group rotation and lattice vibra- 
tions, no other dynamics are present in the considered 
dynamic window, the former contributions are elastic, 
and S(Q,u}) is given by: 



S(Q,u) = e- 2W ^[a coh S(Q)6(u) -f 
(a inc -a^)6(u;)+a^S^(Q,u;)] 



(31) 



with (T co h and <7i nc the total coherent and incoherent 
cross-sections and a^f? the incoherent cross-section of 
the three hydrogens of the methyl group. o- co hS(Q) 
is the coherent static intensity. As mentioned above, 
the Debye- Waller factor e~ 2W ^ gives account for the 
intensity loss with increasing temperature due to vi- 
brations. Since the Debye- Waller factor only enters 
as a scaling factor in the total scattering function, it 
can be removed in the analysis procedure of the ex- 
perimental spectra. Finally, it must be stressed that 
the function S(Q,uj) to be compared with the exper- 
imental spectrum is given by the convolution of the 
total scattering function (31) with the instrumental 
resolution R(Q,uj): 



S(Q,w 



) = f S(Q, 



u>)R(Q, uj — ui')du)' 



(32) 



In the following, comparison between experimental 
spectra and any of the theoretical scattering functions 
given above, will be understood as done after convo- 
lution of the latter with the instrumental resolution. 



IV. METHYL GROUP DYNAMICS IN 
POLYMER SYSTEMS 

A. Rotation Rate Distribution Model (RRDM) 

In this Section we summarize the grounds of the 
RRDM. Though the RRDM was initially developed 
for polymer systems, it must be stressed that all its 
basic assumptions can, in principle, be also applied to 
non-polymeric disordered systems (see Section VII). 
In the framework of the RRDM [25,31-36,88-94], it is 
assumed that the only effect of the structural disor- 
der on methyl group dynamics is to introduce a dis- 
tribution of rotational barriers g(Vz), originated from 
the different local environments felt by the individual 
methyl groups. A non-trivial assumption is that the 
distribution of distances between neighbouring methyl 
groups does not introduce significant coupling forces 
for the smallest distances, which, in particular, would 
modify the single-particle Hamiltonian (9) leading to 
a more complex energy level scheme [95-100]. 
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According to this approach, the spectrum for the 
disordered system is simply obtained as a superpo- 
sition of crystal-like spectra weighted by the distri- 
bution g(V 3 ). Due to the normalization condition 
/ 9{V 3 )dV 3 = I7 the elastic terms in Eq. (31) are 
not affected by the introduction of the distribution. 
Hence, in the framework of the RRDM, S(Q,uj) for 
the glassy system is formally equal to Eq. (31). How- 
ever, the introduction of g(V 3 ) modifies the incoher- 
ent contribution for methyl group dynamics, which is 
given by: 

/>oo 

S%?{Q,u>) = g(V 3 )S^(Q,u;,V 3 )dV 3 , (33) 
Jo 

where each single crystal- like spectrum S^(Q, ui, V 3 ) 
evolves with temperature according to the description 
given in Section III, summarized in Eqs. (14,15,18,23, 
24,25,30). The different parameters controlling this 
evolution are therefore, for each individual methyl 
group: 

i) the tunnelling frequency Ut, first librational energy 
.Eoi, and classical activation energy E A - These three 
quantities are direct functions of the potential barrier 
V 3 , and are obtained by calculating the eigenvalues of 
the Hamiltonian (9). Fig. 7 shows the correspond- 
ing results for the purely threefold case. Highly accu- 
rate numerical relations for the latters, in the range 
V 3 < 3500 K, are given by: 

E A (K) = 0.598K, 105 , (34) 
£ i(meV) = 0.470K, ' 548 , (35) 

hu} t (meV) = 

( V \ 106 
0-655 (^1 + ^J cxp[-(vy4) a5 ]. (36) 

In these equations V 3 is given in units of K. 

ii) the preexponential factors Too, j s , jae and 
lE a Eb f° r the different Arrhenius-like equations 
(15,18,23,24). 

iii) the crossover temperature T C (V 3 ), that determines 
which regime -crossover if T < T C (V 3 ) or classical hop- 
ping if T > T C (V 3 )- governs the dynamics of a given 
methyl group of barrier height V 3 at moderate tem- 
perature T. 



Glassy systems show broad distributions g(V 3 ). For 
that reason, shift, broadening, and merging of the in- 
dividual peaks with the quasielastic component, can- 
not be solved from the spectrum for the glass. On the 
other hand, it is expected that the temperature depen- 
dence of individual spectra is highly sensitive to the 
barrier height. In principle, methyl groups with high 
values of V 3 will reach the classical hopping regime at 
higher T c than those with low values of V 3 . Moreover, 
there is a large set of parameters, namely the preexpo- 
nential factors for shift and broadening, which depend 
on the unknown constants {gu} describing the cou- 
pling between the methyl group and the phonon bath 
[Eqs. (27,28)]. Such preexponential factors cannot 
be easily deduced for a glassy system, in contrast to 
the crystalline case [36,48,50,52-59], where shift and 
broadening of the single peaks can be directly observed 
and analyzed (see Figs. 4, 5). As mentioned above, 
the typical values of observed in crystalline sys- 
tems are in the same order of magnitude (a few meV), 
but 7 S , jae and jE a E b can take a wide range of values. 

All these problems enormously complicate the anal- 
ysis of the spectra in glassy systems. Next it is shown 
that many of the unknown present paramaters, namely 
"fAE, lE a E b and 7 S , controlling the temperature depen- 
dence of the crystal-like spectrum for each individual 
methyl group, can be removed without the need of 
carrying out complex calculations, as it would be re- 
quired by the introduction of a detailed model for the 
coupling of the methyl group to the glass vibrations. 
Instead, a functional relation between the crossover 
temperature and the barrier height, T c — T C (V 3 ), is 
assumed [33-36]. It is found that a consistent de- 
scription of the spectra in all the temperature range is 
possible in terms of only three parameters: the aver- 
age and standard deviation of the barrier distribution 
g(V 3 ), and the preexponential factor for classical 
hopping. 

In the approach introduced by the general version 
of the RRDM [33-36], the crossover temperature T c 
is operationally defined as the temperature where the 
rate for coherent quantum tunnelling becomes equal to 
the rate for incoherent classical hopping, i.e., T(T C ) = 
huJt, or from Eq. (15): 

^ = MToo/fo*)' (37) 

As mentioned above, the classical activation energy 
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Ea and the tunnelling frequency u t are direct func- 
tions of Vz- It has been shown (see Refs. [25,35] 
for a detailed discussion) that can be taken as 
a barrier-independent quantity in a good approxima- 
tion. Therefore, T c defined in this way depends only 
on the barrier height V3 and the preexponential fac- 
tor for classical hopping T^. This latter material- 
dependent parameter can be interpreted as a measure 
of the strength of the coupling of the methyl groups 
to the glass vibrations. For a given methyl group 
of barrier V3, the larger — i.e., the stronger the 
coupling — , the lower will be the T c for the onset of 
classical hopping. Two technical approximations are 
introduced: 

i) For each individual methyl group, the preexponen- 
tial factors jae and jE a E b are taken as equal (in the 
following they will be denoted as 7&). Therefore, 



Tae = ^E a E„ = ■y b cxp(-E 01 /k B T). 



(38) 



ii) The activation energy Eg for the shift of the tun- 
nelling lines is assumed to be equal to E i'. 



KAw t = -y s exp(-E i/k B T). 



(39) 



The general experimental evidence in crystalline sys- 
tems for the values of Jae, lE a E b and Es (see above) 
supports the reliability of these two approximations. 
In addition, two constraints are imposed to Eqs. 
(38,39) at the crossover temperature: 

i) The tunnelling frequency must be shifted to zero 
at the onset of the classical regime, Aw t (7c) = w t , 
as required by the disappearence of quantum effects 
above T c . Hence, from Eq. (39): 



7 S = huj t exp(E 01 /kT c ) 



(40) 



ii) Continuity condition for the HWHM of the 
Lorentzians at the merging point T c , i.e., Tae(T c ) = 
r EaEb (T c ) = r(T c ). Hence, from Eqs. (15,38): 



7b — 1 00 exp[(S i -E A )/kT c }. 



(41) 



A straightforward consequence of Eqs. (37,40,41) is 
that 7,5 =76, which is also in agreement with the gen- 
erally close values found for both quantities in crys- 
talline systems. 



All the quantities at the right sides of Eqs. 
(37,40,41) depend exclusively on V3 and on the barrier- 
independent factor Too. In this way, the preexponen- 
tial factors for the crossover, 76 and j s , can be deter- 
mined unambiguously for each value of V3 in the dis- 
tribution 5(^3). Hence, the temperature dependence 
of the spectra for the glass is modelled exclusively in 
terms of and (7(1/3). Two-parameter (average bar- 
rier and standard deviation) simple distributions are 
used, as Gaussian [25,31-34,36,88-93] functions. In a 
few cases [35,94] two-parameter Gamma functions are 
more adequate. In this way, the RRDM only intro- 
duces three independent parameters for the analysis 
of experimental spectra in all the temperature range. 

In the limit of low and high temperature, all the 
individual methyl groups will perform respectively ro- 
tational tunnelling or classical hopping, and the cor- 
responding crystal-like contributions to the spectrum 
of the glass [Eq. (33)] will be governed respectively by 
Eq. (14) for tunnelling and by Eqs. (15,30) for hop- 
ping. These limits are referred to as the tunnelling and 
hopping limits of the RRDM. At intermediate temper- 
ature the functional relation between T c and V3 intro- 
duced in Eq. (37) allows one to select the correspond- 
ing dynamic regime for each individual methyl group 
— i.e., crossover, governed by Eqs. (25,38,39,40,41), or 
classical hopping, governed by Eqs. (15,30). 

In a first approximation, a Gaussian distribution of 
purely threefold rotational barriers is assumed: 



g(v 3 ) 



1 



exp 



24 



(42) 



with (V3) the average barrier and cry the standard de- 
viation of the distribution. See Fig. 8a for (V3) = 
800 K, ay — 250 K, as an example of parameters for 
which tunnelling and hopping features are observable 
by neutron scattering (see below). 

The functional relations — e.g., Eqs. (34,35,36) — 
between Ea, E^, huj t , and V3, allow for a straight- 
forward transformation between (/(V3) and the corre- 
sponding distributions of activation energies for classi- 
cal hopping, /(Ea), rotational tunnelling frequencies, 
h(hujt), and librational energies, F(Eoi): 



9(V 3 )dV 3 = f(E A )dE A = 
F(E i)dE Qi = -h(hujt)d(huj t ) 



(43) 



In a good approximation, the classical activation 
energy E^ depends linearly on V3 (see Fig. 7c). 
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Hence, the distribution f(E\) is also Gaussian (see 
Fig. 8b), with average energy (£U) and standard devi- 
ation o e- As mentioned above, the preexponcntial fac- 
tor for classical hopping is assumed to be barrier- 
independent. This approximation allows one, accord- 
ing to Eq (15), to transform the Gaussian distribution 
f(E\) into a log-Gaussian distribution of Lorentzian 
HWHM for classical hopping: 



h (log r) 



27TO- 



cxp 



(io g r-io g r ) 



2a 2 



where 



T = Too cM-{E A )/k B T) 



(44) 



(45) 



is the corresponding Lorentzian HWHM for the aver- 
age activation energy, and 



a = a E log(e)/k B T 



(46) 



is the standard deviation for H(\ogT). Fig. 9 shows 
the temperature dependence of this distribution for 
the above parameters (V3) = 800 K, cry — 250 K, and 
for a typical value of Too = 6 meV. As can be seen, at 
high temperature the major part of the distribution is 
accessible in the energy window of neutron scattering 
spectrometers. 

The parameters To and a at a given high temper- 
ature can be obtained from a fitting procedure of the 
experimental methyl group incoherent contribution 
S^(Q, to) to the classical limit of Eq. (33). From the 
transformation g(V3)dVs — — H(logT)d(\og T), and 
from the incoherent scattering function (30) for a sin- 
gle barrier, Eq. (33) can be rewritten in the classical 
limit as: 



1 + 2 jo (Qr) 



2[1 -j (Qr)} 



inc 

/-•DC 



+ 



/CO 
H(logT)L(u;T)d(logT). (47) 
-00 



The set of values To(T) and a(T) obtained for the dif- 
ferent temperatures is fitted to Eqs. (45,46) in order 
to obtain Too, (Ea) and <je- Finally, the knowledge 
of these two latter quantities allows, by transforma- 
tion, to determine (V3) and cry, i.e., the parameters 
of the barrier distribution (7(1/3). The latter, as ob- 
tained from the analysis of classical high temperature 
spectra, can next be transformed into the distribu- 
tion of tunnelling frequencies h(hcj t ) and librational 



energies F(E i) (see Figs. 8c, 8d). The consistency of 
the RRDM requires that these two latter distributions 
derived in this way, also reproduce, respectively, the 
corresponding experimental tunnelling and librational 
spectra. 

Due to the approximate exponential dependence of 
huj t on V3 (see Eq. (36) or Fig. 7a), h{%uj t ) takes 
an extremely asymmetric shape (see Fig. 8d) with 
the maximum shifted to low frequencies. From the 
incoherent scattering function for a single barrier (14), 
the tunnelling limit of Eq. (33) is given by: 



2[l-jo(Qr) 



[h(u t ) + h{-0J t )}. 



(48) 



For sufficiently high average barriers (V3) or broad dis- 
tributions g(Vs), the maximum of h(±hu t ) is placed 
beneath the instrumental resolution (sec Fig. 8d). 
Hence, the tunnelling spectrum shows an apparently 
quasielastic contribution, which could be misinter- 
preted as a signature of hopping events at very low 
temperature. However, its actually inelastic origin 
finds a natural explanation, within the framework of 
the RRDM, in terms of a distribution of rotational 
tunnelling lines. 

Once the tunnelling and hopping limits of the 
RRDM are able to give account, with the same set of 
parameters (V3), cry an d T^, respectively for the ex- 
perimental spectra at very low and high temperature, 
the consistency of the general version of the RRDM 
is tested, with such parameters, at intermediate tem- 
perature. It must be noted that this step does not 
involve any further fitting procedure. Instead, exper- 
imental spectra are directly compared with the the- 
oretical ones, which are constructed according to the 
procedure exposed above for the general version of the 
RRDM, and by making use of the operational defini- 
tion of the crossover temperature, Eq. (37), in order 
to select the dynamic regime (crossover or classical) 
for each individual methyl group. 

Eq. (37) also introduces a distribution of crossover 
temperatures, G(T C ). Fig. 10 shows the latter, for 
the same parameters, (V 3 ) = 800 K, a v = 250 K and 
= 6 meV, of Figs. 8 and 9. As temperature T 
increases, the fraction of methyl groups that reach the 
classical hopping regime, i.e., the area of G(T C ) for 
T c < T, will also increase and will be negligible above 
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some high temperature, where quantum effects will be 
unsolvable from the dominating hopping dynamics. In 
such conditions, the system will be well described by 
the classical limit of the RRDM. 



B. Other approaches 

An alternative approach to the classical limit of the 
RRDM, also based on the idea of an underlying distri- 
bution of rotational barriers, has been introduced by 
Arrighi et al. [23,24] for analyzing methyl group dy- 
namics in the high temperature hopping regime. By 
Fourier transforming into the time domain, the in- 
termediate incoherent scattering function for methyl 
group classical hopping is given by: 



Obviously for a crystalline system F(t) — exp(— t/r), 
with t the single residence time for hopping. 

For a disordered system, the corresponding function 
F(t) in the framework of the RRDM is a log-Gaussian 
distribution of time exponential functions: 



m = f 

J —i 



F(logr)exp(-Vr)d(logT), (50) 



as trivially obtained by Fourier transformation of Eq. 
(47). In the approach introduced by Arrighi et al., 
no functional form is assumed for the distribution 
ff(logr). On the contrary, T{t) is effectively rep- 
resented by a stretched exponential, or Kohlraush- 
Williams- Watts (KWW) function, exp[-(i/T ) /3 ], with 
To a characteristic rotational relaxation time. The 
stretching exponent f3 takes values between and 1. 

A KWW function can be effectively represented as 
a log-distribution of time exponential functions: 



eM-(t/r o ) } 



f 

J —i 



J FC(logr)exp(-t/r)rf(logr). (51) 



The distribution function _ftT(logr), which can be ob- 
tained by inverse Laplace transformation methods 
[101], is asymmetric. Therefore it is clear that the 
KWW function cannot be equivalent to the RRDM 
distribution F(t) if a log-Gaussian form is selected for 
i/(logr) [89], (note that for the same reason tq ^ t). 
This fact is illustrated in Fig. 11 for typical values r 
= 0.1 ns and a = 1. 



Both pictures (KWW and RRDM) achieved descrip- 
tions of similar quality for experimental spectra in a 
limited energy window, i.e., by using a single spec- 
trometer. The corresponding times r and To followed 
Arrhenius behavior with close values of the respective 
activation energies [23,24,89]. However, an analysis in 
a wide dynamic range [89] , by combining several spec- 
trometers, showed that, while an excellent agreement 
of the experimental T[t) with the RRDM function (50) 
was obtained in all the time interval, a simultaneous 
fit to a KWW function was not satisfactory. More- 
over, independent fits in different experimental win- 
dows provided incompatible values of To and (3. This 
result is illustrated in Fig. 12 for the same function 
F(t) of Fig. 11, now analyzed in different narrow dy- 
namic windows. 

Another disadvantage of using a KWW functional 
form for T{t) is that, in contrast to the unified pic- 
ture introduced by the RRDM, it does not provide a 
direct physical connection between the classical and 
quantum features observed respectively at high and 
low temperature. It is worthy of remark that MDS 
support, by following separately the behavior of each 
individual methyl group in the system, the general va- 
lidity of the Gaussian approach for the distribution of 
rotational barriers [85,86,88,102]. 



V. APPLICATION OF THE RRDM TO A 
SHOWCASE: POLY(VINYL ACETATE) 

PVAc was the first polymer — and indeed the first 
highly disordered system — where a complete inves- 
tigation by neutron scattering was carried out on 
all the features of methyl group rotation (tunnelling, 
crossover, hopping and librations). It was also the first 
system where a fully consistent description of the ex- 
perimental results was achieved in terms of the RRDM 
[31,33,89]. Figs. 13, 14 show respectively for a BS and 
a TOF spectrometer, experimental spectra (circles) for 
PVAc, at several temperatures in the tunnelling (T = 
2 K), crossover (T = 20 and 30 K) and classical hop- 
ping (T = 70, 120, 160 and 200 K) regimes. All these 
temperatures are far below the glass transition tem- 
perature (T g « 315 K), where motions different from 
methyl group rotation just contribute as vibrational 
dynamics via the Debye- Waller factor. 

An excellent agreement with the RRDM is observed 
in all the temperature range (solid lines). As in the 
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usual analysis procedure, high temperature data of 
PVAc were first analyzed in terms of the classical 
RRDM [89]. Fig. 15 shows the temperature depen- 
dence of the parameters log To and a of the distribu- 
tion H(\ogT) (taken as Gaussian) obtained from the 
analysis in terms of the hopping limit of the RRDM. 
The use of several spectrometers with different energy 
resolutions and dynamic ranges reduced the uncertain- 
ties of the so determined parameters. The observed de- 
viations below T w 70 K from the linear behavior pre- 
dicted by Eqs. (45,46) show the relevance of quantum 
effects at low temperatures. The fact that the mag- 
nitude of such deviations depend on the experimental 
dynamic window evidences the importance of using 
several spectrometers in order to avoid, below some 
temperature, an analysis biased by the used instru- 
mental resolution. The results of the fits for T > 70 
K to Eqs. (45,46) provide a preexponential factor 
Toe =9.1 meV, an average classical activation energy 
(Ea) = 450 K, and a standard deviation oe — 250 K. 
Transformation to g{V^) yields (V3) = 534 K and oy = 
274 K. The dashed lines in Fig. 15, corresponding 
to r = Foo cxp[— {{Ea) ± cte)/^bT], have been intro- 
duced in order to stress the presence of a wide distri- 
bution of hopping times, which progressively spreads 
over several orders of magnitude as temperature de- 
creases. 

The corresponding distribution of tunnelling fre- 
quencies h{Tuot) is obtained by transformation of the 
barrier distribution g{V^) independently derived from 
the high temperature classical analysis. Once h(haj t ) 
is known, the theoretical S^(Q, u) is constructed ac- 
cording to Eq. (48). An excellent agreement between 
the tunnelling limit of the RRDM and the experimen- 
tal spectrum at 2 K is obtained (Fig. 13). The max- 
imum of h(%ujt) is at hu>t ~ 0.03 /zeV [31], i.e., well 
beneath the instrumental resolution, leading to the ap- 
parently quasielastic observed intensity . 

For the case of the librational distributions, a direct 
comparison with experimental results is complicated, 
since librational energies are strongly superposed with 
the phonon spectrum. This is particularly problem- 
atic in PVAc, which cannot be selectively deuterated 
in order to attenuate the intensity scattered by nuclei 
different from methyl group protons. Fig. 16 shows a 
comparison of the generalized VDOS of PVAc with the 
distribution of librational energies F(E i) obtained by 
transformation of g(Vs). A reasonable agreement be- 



tween the maximum of the experimental peak and 
the maximum of F(Eo\) is achieved, supporting the 
validity of the previuosly obtained distribution g(Vs) 
and the consistency of the RRDM picture. It must 
be stressed that, since the generalized VDOS is not 
the real density of states (see Section II), a rigorous 
comparison between the widths of the experimental 
librational peaks and the theoretical librational distri- 
butions cannot be made. 

Once it has been checked that the distribution of 
potential barriers g(Vs) can give account, by transfor- 
mation to the other distributions in Eqs. (43, 44), 
for the spectra in the tunnelling and hopping dynamic 
limits, as well as for the librational peak observed in 
the generalized VDOS, the consistency of the general 
RRDM is confirmed by also reproducing the experi- 
mental spectra at intermediate temperatures (Fig. 13 
for T = 20 and 30 K). An excellent agreement is again 
achieved. A more detailed comparison between exper- 
iment and theory can be obtained by investigating the 
temperature dependence of the integrated intensities 
in different inelastic windows, as shown in Fig. 17. 
Since Eq. (48) is temperature-independent, and the 
Dcbye- Waller factor is a decreasing function of T , the 
initial increase of the integrated intensities when heat- 
ing the system makes clear that the tunnelling limit 
of the RRDM is not appliable for T > 2 K. It must 
be stressed that the observed double-peak structure 
is not an experimental artifact. Indeed, it can also 
be observed in a crystalline system when the inelas- 
tic window is properly selected, e.g., at lower energy 
but close to the single tunnelling line. As mentioned 
in Section III, when increasing temperature the tun- 
nelling line broadens and shifts to lower energies un- 
til it merges, when reaching the hopping regime, into 
a single quasielastic line. As a consequence, the in- 
tensity measured in the selected inelastic window will 
pass through a first maximum. A second maximum 
will be observed at higher temperature when broad- 
ening of the quasielastic line reaches the selected in- 
elastic window. In glassy PVAc, due to the broad dis- 
tribution of tunnelling frequencies in comparison with 
the energy window accesible by the spectrometer, the 
double-peak structure is observed at any inelastic win- 
dow. 

As shown in Fig. 17, the double-peak structure 
is nicely reproduced by the theoretical curves (solid 
lines) obtained by integration of the RRDM scattering 
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function in the corresponding inelastic windows. The 
theoretical curves quantitatively reproduce the exper- 
imental intensities in all the temperature range. It 
must be noted that the formers have been modulated 
by a Debye-Waller factor e - 2 Q 2 <" 2 >/ 3 , with (u 2 ) = 9T 
and = 3 x 1(T 4 A 2 !^ 1 [31]. The theoretical curves 
corresponding to the classical hopping limit of the 
RRDM (dashed lines) have been extrapolated to very 
low temperature in order to stress the relevance of 
quantum effects below T w 70 K, as was pointed out 
above. 



VI. SUMMARY OF EXPERIMENTAL 
RESULTS BY NEUTRON SCATTERING ON 
POLYMERS 

A. Classical hopping 

From early 80 's, neutron scattering has been widely 
used to investigate methyl group dynamics in the 
quasielastic energy range \hu\ < 2 meV, or by us- 
ing the FEW technique (see Section II). A huge col- 
lection of data has been obtained for several poly- 
mer systems, including standard materials as PVAc 
[89,103], PMMA [20,22-24,34,91,104], polypropylene 
oxide) (PPO) [105], polyvinyl methyl ether) (PVME) 
[23,25,27,106], 

polyisoprcnc (PI) [26,92], neat poly(dimcthyl siloxane) 
(PDMS) [93,107-109], and filled with silica nanopar- 
ticles [110], head-to-tail [88,90], and head-to-hcad 
(hh) [111] polypropylene (PP), poly (ethylene propy- 
lene) (PEP) [111], and poly (methyl phenyl siloxane) 
(PMPS) confined in silica based nanoporous glasses 
[112], polymer blends as solution chlorinated polyethy- 
lene (SCPE)/PMMA [22,113], polystyrene (PS)/ 
PVME [114,115], hh-PP/PEP [111], or polyethylene 
oxide) (PEO)/PMMA [116,117], thin film photoresist 
polymers [118], and other polymer based compounds 
as polymer electrolites [119-121], proteins [122], liq- 
uid crystalline polyesters and polyethers [123,124] or 
propylene glycol oligomers in the bulk state [125] and 
confined in clay [126,127]. Very recently, Russina et 
al. have investigated methyl group classical hopping 
in a neat, and doped with an inorganic salt, PEO- 
PPO random copolymer [128]. Differently from the 
usual FEW technique, data have been acquired and 
analyzed in an inelastic fixed energy window centered 
at huj w -15/xeV with BS resolution. 



A simple way to investigate the temperature depen- 
dence of methyl group rotation is by following the tem- 
perature evolution of the FEW intensity, Ifew{Q,T) 
(see Section II). If the corresponding energy window 
for data acquisition is — A e ; < hui < A e ;, then integra- 
tion of Eq. (47), together with the additional clastic 
contributions (see Section III), provides the RRDM 
function for the FEW intensity: 



Ifew(Q,T) = e- 2W( Q- T \<r coh S(Q) + 

(Q,T)), 
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where the incoherent contribution for methyl group 
classical hopping is given by [93]: 

rMG (n ji\ _ 1 + 2j (Qr) 2[l-j (Qr)] 
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J — ( 



H(Iogr)^arctg^) (fog] 



If temperature is sufficiently low so that most of 
the methyl groups rotate with a hopping rate T <C 
A el , the limit (2/7r)arctg(A e; /r) — > 1 is fulfilled and 
Ip^EWiinc^QiT) = 1. In such conditions the temper- 
ature dependence of the FEW intensity is just con- 
trolled by the Debye-Waller factor, Ifew(Q,T) = 
e -2W(Q,T) [ CTcoh 5'(Qj + CTinc ]. it must be reminded that 

below some temperature the classical hopping picture 
is not valid, and quantum effects dominate methyl 
group dynamics. However, from integration of Eqs. 
(48) and (25,33) in the FEW window, the latter result 
can be extrapolated to very low temperature if, for 
most of the methyl groups, Tae, r ' E a E b and Vt are & U 
smaller than A e ;. 

FEW data for several polymers [26,90,93], obtained 
by means of high resolution BS spectrometers, are 
shown in Fig. 18. Intensities are normalized to their 
values at T — > 0. The onset of methyl group dynamics 
with time scales shorter than that corresponding to 
the instrumental resolution (~ 5 ns for the used spec- 
trometers) manifests as a step- like decrease of the nor- 
malized intensity, which, above some temperature, de- 
viates from the exponential (linear in the logarithmic 
scale of Fig. 18) decrease corresponding to the Debye- 
Waller factor, i.e., to pure vibrational dynamics. The 
second step at the very high temperature range corre- 
sponds to the onset of much slower processes such as 
secondary relaxations or motions involved in the glass 
transition process. 
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The FEW technique requires much less acquisi- 
tion time than QENS measurements, but obviuosly 
quasielastic spectra are preferred for an accurate de- 
termination of the parameters of the distribution of 
hopping rates. Though for moderate average barri- 
ers, analysis of FEW data have provided reliable re- 
sults [90], this procedure is not recommended when 
the average rotational barrier is very high (V3 > 2000 
K), since in such cases the two step- like decays of the 
FEW intensity overlap. Such a case is illustrated for 
polysulfone (PSF) in Figure 18. 

First quasielastic investigations on methyl group 
classical rotation in polymers were analyzed in terms 
of the usual approach for molecular crystals, i.e., by us- 
ing a single Lorentzian function (single hopping rate) 
to model the quasielastic component. As mentioned in 
the Introduction, such a procedure fails when applied 
to polymers (and in general to highly disordered sy- 
sems), providing inconsistent and unphysical results. 
In polymeric materials, and in general in disordered 
systems [35,36], the effect of structural disorder has 
to be taken into account to describe properly methyl 
group dynamics. In the framework of the RRDM, it 
is assumed that the only effect of the structural disor- 
der on methyl group dynamics is to introduce a dis- 
tribution of rotational barriers giV^), originating from 
the different local environments felt by the individ- 
ual methyl groups. The initial version of the RRDM 
— the hopping limit — was introduced for analyzing 
quasielastic data from P VME [25] . Nearly at the same 
time, a very close approach based in a Gaussian dis- 
tribution of activation energies was proposed by Frick 
and Fetters [26] to describe the temperature depen- 
dence of FEW data from PI. As mentioned in Section 
IV, an alternative approach [23,24], first introduced 
to analyze time-domain data of PMMA and PVME, 
assumed a KWW stretched exponential function for 
the intermediate scattering function I(Q,t). Although 
this procedure was able to describe spectra obtained 
in a single spectrometer, it failed to properly account 
for data collected in a wide dynamic range by using 
several instruments [89] . 

Classical hopping between three equivalent posi- 
tions was assumed in all these approaches. It must be 
noted that classical hopping is just driven by the acti- 
vation energy Ea, and that a purely threefold barrier 
and another one distorted by a small sixfold contribu- 
tion, T4(l — cos(60 + 6))/2, can provide the same Ea 



with an appropiate selection of the parameters V3, Vq 
and 5. Hence, an analysis of high temperature classi- 
cal spectra provides a first approach to the functional 
form of the rotational potential, but an unambiguous 
determination of the latter is only possible by means of 
a complete analysis of hopping, tunnelling and libra- 
tional spectra. For most of the polymers or molecular 
glasses for which at least two of these three features 
have been investigated, it has been found that a dis- 
tribution of purely threefold barriers provides a good 
description of experimental spectra. As shown below, 
small but significant corrections to the threefold term 
have instead to be considered in order to consistently 
reproduce the different features of the ester-methyl 
group dynamics in PMMA. 

Recent QENS investigations in PDMS [93], PVAc 
[103] and in the blends hh-PP/PEP [111], and 
PEO/PMMA [116] have tackled the problem of methyl 
group classical hopping in the very high temperature 
range where secondary relaxations or motions involved 
in the glass transition process enter the dynamic win- 
dow of neutron scattering. In that situation the time 
scales of such processes and methyl group hopping su- 
perpose. A consistent description has been achieved 
by assuming that the mentioned motions and methyl 
group rotation are statistically independent. This ap- 
proach was first introduced in an investigation, by 
means of molecular dynamics simulations, on methyl 
groups dynamics in PI above the glass transition tem- 
perature [129]. Within this approach, the RRDM pic- 
ture is maintained for methyl group rotation and the 
total scattering function is just obtained as a convo- 
lution of the RRDM scattering function (47) and that 
corresponding to the other relaxational modes [129]. 
A unified description is achieved with the same values 
— or small variations — of the RRDM parameters de- 
rived at lower temperatures, where the dynamic con- 
tribution of the polymer matrix just enters in the vi- 
brational Debye- Waller factor. This result evidences 
that the barrier distribution for methyl group rotation 
in polymers is not substantially changed by the struc- 
tural rearrangements produced when approaching the 
glass transition [93,129]. 
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B. Rotational tunnelling and crossover to hopping 

As mentioned in the Introduction, the existence of 
quantum rotational tunnelling for methyl group dy- 
namics in polymer systems was expected for a long 
time but it was not directly observed by neutron scat- 
tering until late 90's [31,32]. In contrast to crystalline 
systems where, due to the unique value of the ro- 
tational barrier, resolution-width tunnelling lines are 
observed at very low temperature, in amorphous sys- 
tems, a distribution of rotational barriers, and there- 
fore of tunnelling frequencies, is expected. Some ex- 
amples are misoriented methyl groups in crystalline 
lattices [130] or absorbed in water clathrates [131] or 
crystalline zeolites [132], and other small molecular 
rotors in chemically or geometrically disordered en- 
vironments as non-stochiomctric metal hexaammines 
[133], ammonia in metal alkali fullerides [134], diluted 
solutions of methane in noble gases [135-137], or of 
ammonium in metal alkali halides [138,139]. In these 
systems, disorder effects provide narrow distributions 
of rotational barriers and of the corresponding tun- 
nelling frequencies, which just manifest as weak broad- 
ening of the tunnelling peaks in neutron scattering 
spectra. However, in the case of methyl groups in 
polymers [31,32,34], molecular glasses [35,36], or close 
to pore walls in mesoporous silicates [94,140-142], dis- 
order leads to broad distributions of tunnelling fre- 
quencies. Due to the typical values of the rotational 
barriers in the latter systems, specially in polymers, 
most of the expected distribution of tunnelling lines 
corresponds to energies below ~ 1 /xeV, i.e., beyond 
the dynamic limit explored by BS spectrometers. This 
fact prevented for a long time the observation of tun- 
nelling features in low temperature spectra. 

PVAc is the polymer having the lowest energy bar- 
rier for methyl group rotation reported so far, and the 
most suitable one for the direct observation of rota- 
tional tunnelling transitions by means of neutron scat- 
tering [31,33]. An analogous investigation has been 
carried out on rotational tunnelling of the ester-methyl 
group of PMMA [32,34]. Though having a higher bar- 
rier, PMMA has the advantage of being selectively 
deuterablc, in contrast to PVAc. In this way, the ob- 
served scattered intensity excess for PMMA-d5 (i.e., 
with partial deuteration of the a-methyl group and 
the main-chain hydrogens) over the elastic line at very 
low temperature can be unambiguously assigned to the 



ester- methyl group dynamics. 

A further support of the interpretation of the ap- 
parently quasielastic low-temperature intensity as the 
result of a distribution of rotational tunnelling lines 
is provided by exploting the well-known isotope ef- 
fect. Hence, deuteration of the methyl group yields 
a rotational constant B twice smaller than for the 
protonated methyl group. The corresponding tun- 
nelling frequency is reduced by several orders of mag- 
nitude (dashed line in Fig. 7a), and consequently, the 
corresponding distribution h{TitOt) is strongly shifted 
to lower frequencies (see Fig. 19 for the ester-methyl 
group of PMMA). This effect predicts a strong su- 
pression of the scattered intensity excess in the exper- 
imental spectrum, as it is indeed observed [32] in a 
fully deuterated sample (d8) of PMMA at the same 
temperature (see Fig. 20). For a proper compari- 
son between the ester-methyl group dynamics in the 
d5- and d8-samples, the elastic contributions, corre- 
sponding to the total coherent cross-section and to 
the incoherent cross-section of the other nuclei, are 
substracted from the spectra. This procedure intro- 
duces large uncertainties in the so obtained spectrum 
for the d8-sample, due to the large relative weight of 
such elastic contributions. However, the supression in 
the d8-sample of the apparently quasielastic intensity 
is evident, and consistent with the result predicted on 
the basis of the RRDM parameters derived for the d5- 
sample [32]. 

Except for PVAc and for the ester-methyl group of 
PMMA, for the rest of the polymers investigated so 
far, the large values of the average barriers, (V3) > 
1000 K, determined from the analysis of high tempera- 
ture hopping data and/or of librational levels, prevent 
the observation of clear rotational tunnelling features 
at very low temperature. This situation is illustrated 
in Fig. 21 by a comparison between the expected dis- 
tributions of tunnelling frequencies for several poly- 
mers. As can be seen, within the neutron scattering 
window, h(hu>t) for PVME and PI is orders of magni- 
tude smaller than for PVAc and for the ester-methyl 
group of PMMA. Convolution of h(Uui t ) with the in- 
strumental resolution provides a scattering intensity 
excess unsolvable from the latter within the experi- 
mental noise, as indeed is experimentally confirmed. 
An analysis of high temperature data for PDMS in 
terms of the hopping limit of the RRDM [93] provided 
a distribution of classical activation energies J(Ea), 
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which, by transformation, would yield a distibution of 
tunnelling frequencies clearly observable by means of 
high-resolution BS spectrometers. However, the mea- 
sured intensity at T — 2 K [143] is hardly distinguish- 
able from the instrumental resolution, as can be seen 
in Fig. 22. Moreover, the RRDM parameters given 
in Ref. [93] provide, by transformation, a theoretical 
librational maximum around sa 17 meV, far from the 
experimental value of 21.2 meV (see Fig. 23a). On 
the contrary, the RRDM parameters in Ref. [89] re- 
produce the librational peak and are compatible with 
the experimental spectrum at T = 2 K (see Fig. 22). 
The reason for the mentioned uncompatibilities with 
parameters in Ref. [93] remains to be understood. 

Despite of the existence of very few polymer sys- 
tems where methyl group rotational tunnelling can be 
observed by neutron scattering, successful investiga- 
tions have been carried out in two molecular glasses 
— toluene [35] and sodium acetate trihydrate [36]— 
supporting the validity of the RRDM picture for rota- 
tional tunnelling. 

As mentioned in Section IV, at the crossover inter- 
mediate temperature regime (typically 20 K < T < 
100 K) quantum effects arc relevant for methyl group 
rotation, and quasielastic spectra cannot be repro- 
duced in terms of the classical limit of the RRDM. 
However, the fact that spectra do not evolve with tem- 
perature according to the simple expectation for vibra- 
tional dynamics (see Section V) also evidences that 
the tunnelling limit of the RRDM cannot cither be 
applied, and making use of the general version of the 
model is neccessary. The crossover regime has been 
investigated in detail in the only two polymers (PVAc 
[33] and PMMA [34]) where the distribution of rota- 
tional tunnelling lines can be (partially) explored by 
neutron scattering. The success of the general version 
of the RRDM for reproducing all the temperature evo- 
lution of the spectra for both polymers has been ex- 
tended to the glassy states of toluene [35] and sodium 
acetate trihydrate [36]. 

C. Librational levels 

First neutron scattering experiments on methyl 
group dynamics were focused on the determination 
of the librational levels. Since early 70's INS inves- 
tigations on methyl group librations have been per- 



formed in a huge collection of polymeric materials as 
PPO [16], PP [144,145], PMMA [16,17,20,34], poly(a- 
methylstyrene) (PaMS) [17], PVAc [146], PDMS [146], 
PI [26,147], PIB [145,147,148], bisphenol-A based en- 
gineering thermoplastics as PSF, polycarbonate (PC) 
and phenoxy (PH) [149], and layered polytypes of hex- 
atriacontane [150]. 

Figs. 16, 23, 24 and 25 show some examples of 
the experimental generalized VDOS for PVAc [146], 
PDMS [146], PMMA (ester-methyl group) [34], PI 
[26] and PSF [149], compared with the theoretical 
librational distributions F(Eq\) obtained within the 
RRDM approach. For moderate and low energy bar- 
riers (V3 < 500 K), the lowest librational transitions 
occur in the same energy range (5 — 15 meV) where 
a broad shoulder in the generalized VDOS is observed 
for any disordered system [151]. This strong super- 
position complicates the identification and analysis of 
the corresponding librational peaks, as in the case of 
PVAc [146] (Fig. 16). On the contrary, for moder- 
ate and high rotational barriers (V3 > 1000 K) as in 
PDMS, PI or PSF, librational peaks are well resolved 
from the generalized VDOS. 

In order to solve methyl group librations, it is use- 
ful to carry out and compare INS measurements in 
samples with both protonated and deuterated methyl 
groups (see Figs. 24 and 25 for PI and PSF respec- 
tively). It must be stressed that even in the case of 
partial deuteration of all the hydrogens not belonging 
to the methyl group, identification of librational peaks 
is still a complicated task, since any other vibrational 
mode involving the motion of the methyl group hy- 
drogens will provide a peak in the INS spectrum. A 
comparison with infrared and Raman measurements 
is helpful for a differcntation of the physical origin of 
the corresponding peaks. A nice example has been re- 
cently reported by Adams et al. for PIB and PI [147]. 

A relevant point for the analysis of librational fea- 
tures in INS spectra is the high sensitivity of the li- 
brational energies to the shape of the rotational bar- 
rier. Hence, the introduction of moderate sixfold cor- 
rections to the main threefold term of the rotational 
potential provides a significant shift of the librational 
energy. Figure 23b shows how the introduction of a 
small sixfold contribution is sufficient to account for 
the first librational peak of the ester-methyl group in 
PMMA [34], by using the same RRDM parameters 
describing both the rotational tunnelling and classical 



22 



hopping regimes. It is noteworthy that such a con- 
sistent description has been achieved by fixing a ra- 
tio V0/V3 = 0.11 for all the individual methyl groups. 
Hence, as in the purely threefold case, the only inde- 
pendent parameter affected by the distribution effects 
is the threefold term V3. 

Librational energies are also sensitive to coupling 
between pairs of neighbouring methyl groups. In such 
cases, the single-particle approach is modified by in- 
troducing a coupling Hamiltonian. In the simplest ap- 
proximation the coupling term is purely threefold and 
no angular offsets are included: 



& 2 



d 2 



+y(l-cos3^) + -^-(l-cos3(^-V)). (54) 

In a recent INS investigation on methyl group libra- 
tions in PSF [149], the double-peak structure in the 
range w 34 — 44 meV (see Fig. 25) has been inter- 
preted as the result of a distribution of splitted first 
librational energies, such a splitting resulting from 
the coupling between the two methyl groups in each 
monomeric unit. As in the single-particle approach, 
the terms U3 and V3 are taken as Gaussian distributed. 
Due to the highly disordered character of the poly- 
mer matrix and the symmetric chemical arrangement 
of both methyl groups in the monomeric unit, it is 
expected that, on average, they feel the same typi- 
cal environment. Hence, the distributions g(Uz) and 
g(Vs) are assumed to be identical. Since the coupling 
term W3 is strictly due to the interaction between the 
two methyl groups of the same monomeric unit, it is 
expected to be weakly affected by local packing con- 
ditions. Hence for simplicity, distribution effects are 
neglected for the coupling term. 

The classical hopping regime for PSF is not accesi- 
ble by QENS at low and moderate temperature, and 
at high temperature superposes with secondary relax- 
ations which complicate the analysis. Instead, it can 
be properly explored by NMR techniques. Hence, the 
temperature dependence of the D-NMR line shape of 
PSF has been reproduced in terms of a log-Gaussian 
distribution of classical hopping frequencies as Eq. 
(44), with logTo and a respectively following Eqs. 
(45, 46). Analogously to the procedure exposed in 
Section IV, these results provide a Gaussian distri- 
bution of classical activation energies J(Ea), which 



can be transformed into a distribution of rotational 
barriers g(Vz) and librational energies F(E i) via the 
functional relations Ea = £^4(^3), E i = £01(^3). 
Such functional relations can be obtained by solving 
the eigenvalues of the Hamiltonian (54) for the selected 
value of W§ and a large set of values U3 = V3 (sec Rcf. 
[149] for numerical relations). 

Fig. 25 shows the corresponding distributions 
F(Eqi) for the single-particle approach, and for cou- 
pling of methyl groups, — with the conditions above 
explained for the second case. While the single- 
particle approach only gives account for the lowest ex- 
perimental peak, the introduction of a coupling term 
reproduces the double-peak structure. Analogous re- 
sults are obtained for PC and PH [149]. It is note- 
worthy that the introduction of a coupling term does 
not substantially complicate the original RRDM, since 
the only model parameter affected by the distribution 
effects is, as in the single-particle case, the rotational 
barrier height V3. 



VII. MOLECULAR ORIGIN OF THE 
BARRIER DISTRIBUTION FOR METHYL 
GROUP DYNAMICS 

The results summarized in Table II show that the 
average barrier height (V3) for methyl group reorien- 
tation in amorphous polymers is highly sensitive to 
the chemical structure of the monomeric unit. This 
fact suggests that the rotational potential is only par- 
tially determined by intcrmolccular interactions. On 
the other hand, there is not a clear correlation between 
the width of the barrier distribution and the chemical 
structure of the monomeric unit. Next we summarize 
the most relevant neutron scattering results relating 
the parameters of the barrier distribution with pack- 
ing conditions and intermolecular interactions. 



A. Influence of the chain conformation 

The relative orientation between identical side 
groups placed at different monomeric units can be 
random or can adopt specific forms. Hence, in the 
isotactic (i) and sindiotactic (s) conformations, side 
groups at adjacent monomeric units are respectively 
paralell and anti-paralell. If both relative orientations 
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are randomly distributed, the conformation is hetero- 
tactic (h). When there is no dominating tacticity in 
the ensemble of chains forming the system, the latter 
is atactic (a). The ester- methyl group dynamics in 
PMMA have been analyzed in terms of the RRDM for 
different degrees of tacticity of the ester group. The 
influence of local packing conditions on the barrier dis- 
tribution for samples chemically identical is evidenced 
by the small, but significant, differences between the 
obtained RRDM parameters [32,34,91,152] (see Ta- 
ble II). Hence, the average activation energy for a 
purely syndiotactic sample is (Ea) = 710 K [34], while 
for an atactic sample with relation s:i:h — 50:10:40, 
(E A ) = 529 K [91]. 

PMMA chains can be rearranged to form a stereo- 
complex form, where syndiotactic chains wrap around 
isotactic ones to form a double stranded helix struc- 
ture. In Ref. [91] stereocomplexed samples with com- 
plementary deuteration of chains with different tac- 
ticities were investigated. In this way, the effect of 
stereocomplcxation on methyl group dynamics in dif- 
ferent types of chains could be discriminated. For the 
isotactic chains, stereocomplexation considerably re- 
duced the width of /(Ea), yielding a standard devia- 
tion <te = 205 K, as compared to the value oe = 313 
K for the non-stereocomplexed sample [91]. 

For the case of the a-methyl group in PMMA, INS 
measurements showed a much higher librational peak 
[17], and consequently a much higher average rota- 
tional barrier, for the syndiotatic form as compared 
with the isotactic one. On the contrary, INS spectra 
for PaMS [17] did not show significant differences be- 
tween a syndiotactic sample and an atactic one. FEW 
scans reported in Ref. [153] for different tacticities of 
PP did not show apparent differences, within the ex- 
perimental noise, in the temperature range dominated 
by methyl group dynamics. This result is supported 
by the rather close values of the measured librational 
peaks for the isotactic [144] and atactic [145] confor- 
mations. 

There are very few systems where experimental data 
are available for both, the usually investigated head- 
to-tail, and the head-to- head sequences of adjacent 
monomeric units. For PaMS [17] and for PP [111] 
it has been found, repectively by means of INS and 
QENS measurements, that the average barrier in the 
head-to-head conformation is much lower than in the 
head-to-tail one. In the case of hh-PP, the distribution 



of classical activation energies /(Ea) has been calcu- 
lated from an analysis of high temperature data in 
terms of the hopping limit of the RRDM [111]. Appart 
form the mentioned decrease of the average barrier, 
the analysis provides a much narrower distribution 
than that obtained for head-to-tail PP [88,90]. It is 
worthy of remark that the maximum of the librational 
distribution F(E i) for hh-PP, obtained by transfor- 
mation from /(Ea) in a purely threefold approach (see 
Section IV), is close to the experimental value reported 
in Ref. [145], though the agreement is not fully satis- 
factory [111]. A possible origin of such a difference 
might be, similarly to the case of the ester-methyl 
group in PMMA (see Section VI), the presence of small 
higher-order corrections to the main threefold term of 
the rotational potential, or small coupling effects be- 
tween neighbouring methyl groups, which could shift 
the librational energy to slightly higher values than 
the expected ones for a purely threefold single-particle 
potential. This hypothesis is consistent with the prox- 
imity between methyl groups of adjacent units in the 
hcad-to-head conformation of PP. 

B. Influence of mixing with other materials 

A few works report comparisons between methyl 
group rotation in the neat polymer and blended with 
other polymers. QENS investigations in the blends 
PVME/PS [114,115] and hh-PP/PEO [111] did not 
show, for similar concentrations of the two compo- 
nents, changes within the error bar on the barrier dis- 
tribution as compared respectively to neat PVME and 
hh-PP. High dilution of PVME in PVME/PS provided 
a significantly broader distribution, but again did not 
affect to the average barrier [115]. On the contrary, 
progressive dilution of PMMA in SCPE/PMMA in- 
creased the average barrier for the ester-methyl group 
even at moderate concentrations of SCPE [114]. 

The effect of doping with small molecules on methyl 
group dynamics has also been reported in the litera- 
ture, though a detailed quantitative analysis in terms 
of barrier distributions is lacking. Hence, FEW scans 
in Ref. [110] showed no significant differences for neat 
PDMS and filled with silica nanoparticles. Some inves- 
tigations on polymer electrolites suggest that methyl 
group dynamics are not significantly affected by addi- 
tion of the salt [119,120,128], but the present informa- 
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tion is rather limited to make reliable conclusions at a 
quantitative level. 

C. Molecular glasses and comparison with the 
crystalline state 

It is worthy of remark that, not only amorphous 
polymers, but also the low-molecular- weight glasses 
investigated until now display a similar phenomenol- 
ogy for methyl group dynamics. In contrast to amor- 
phous polymers, glass-forming molecular systems can 
be easily obtained in the crystalline state. This possi- 
bility allows one to make a direct comparison between 
the methyl group dynamic features in the crystalline 
and glassy states of a same system. This comparison 
was done in situ for a same sample of sodium acetate 
trihydrate [36]. The analysis of the neutron scatter- 
ing spectra, in the crystalline and glassy states, unam- 
biguously demonstrated that the apparently quasielas- 
tic intensity observed for methyl group tunnelling in 
the glassy state is a consequence of the atomic dis- 
order. Fig. 26 shows the corresponding tunnelling 
spectra for the crystalline and glassy state of sodium 
acetate trihydrate. For this system it was found that 
the average potential barrier in the glassy state takes, 
within the experimental error, the same value as the 
unique barrier in the crystal. A direct comparison be- 
tween methyl group dynamics in the crystalline and 
glassy state has also been reported for ethylbenzene 
[154]. Disorder effects are evidenced by the broad 
structure of the librational peak in the glassy state 
as compared to the resolution-width line in the crys- 
talline case. The small shift of the experimental max- 
imum indicates that, as in sodium acetate trihydrate, 
the distribution is centered around the unique barrier 
of the crystal. On the contrary, an investigation on 
glassy toluene has reported an average rotational bar- 
rier much larger than the unique one of the crystalline 
state [35]. 

The widths of the barrier distributions obtained for 
the former molecular glasses take values comparable to 
those obtained for polymer glasses. All these results 
strongly suggest that the origin of the distribution of 
potential barriers in polymers is directly related to the 
structural disorder and packing conditions, and not 
with chain connectivity aspects which are specific of 
polymer glasses. This hypothesis has been recently 



supported by measurements of the tunnelling spec- 
trum of tri(vinyl acetate) [143]. As can be seen in 
Fig. 27, the experimental spectrum at T — 2 K is re- 
produced by using the same RRDM parameters as in 
PVAc, suggesting that the barrier distribution is not 
affected by polimerization of the basic molecular unit. 

D. Molecular dynamics simulations 

A deep understanding of the microscopic origin of 
the distribution of potential barriers is not easily in- 
ferred from experiments. However, fully atomistic 
MDS arc a suitable tool to gain insight into this prob- 
lem. In particular, the calculation of the distribu- 
tion of librational energies is straightforward and does 
not require lengthy simulation runs. As has been ex- 
posed in previous Sections, once this distribution is 
obtained, the corresponding distribution of potential 
barriers can be directly derived by transformation of 
the former. In addition, MDS offer the possibility of 
varying the parameters of the employed force field, 
and to investigate the influence of the modified force 
field on methyl group dynamics. A comparison be- 
tween MDS and neutron scattering results on PI was 
reported in Ref. [87]. From the atomic trajectories 
computed in the simulations, the generalized VDOS, 
Z(u>), was calculated as the spectral density of the ve- 
locity autocorrelation function: 
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where brackets denote ensemble average over all the 
atoms. The weighting factor g contains, for each 
atomic species, the corresponding cross-section and 
mass, in order to make a direct comparison with the 
experimental results. Fig. 28 shows a comparison 
between the INS data reported by Frick and Fetters 
[26] for PI-d5 (deuterated main-chain and protonated 
methyl group) and the generalized VDOS, computed 
from the MDS as in Eq. (55). A good agreement is 
obtained at least in the energy range, Tiu < 40 meV, 
experimentally investigated in Ref. [26] , validating the 
employed force field. Although quantum effects (not 
considered in the classical MDS) should certainly af- 
fect the high energy range of Z(u>), it is worthy of 
remark that the energies of the computed high energy 
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maxima (see Fig. 28) roughly correspond to the re- 
ported infrarred and Raman bands [155], as well as to 
recent INS results by Adams et al [147]. 

It must be reminded that, although the generalized 
VDOS for PI-d5 is dominated by motions of methyl 
group protons, it contains not only the contributions 
of the librational modes, but also from any other mode 
involving the motion of the methyl group. This is evi- 
denced in Fig. 24, where the librational peak is super- 
imposed to a broad low-energy hump similar to that 
present in the generalized VDOS of PI-d3 (deuterated 
methyl group and protonated main-chain). However, 
the methyl group librational contribution can be sepa- 
rated from the computed generalized VDOS by calcu- 
lating the spectral density of the angular correlation 
function [87]: 
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where brackets denote ensemble average over methyl 
group protons. F(E i) is formally the librational den- 
sity of states. Finally, the distribution of potential bar- 
riers g(Vs) is obtained by transformation from F(Eqi) 
as exposed in Section IV. 

In order to investigate the relation between the aver- 
age and width of the barrier distribution and the non- 
bond interactions, a series of MDS were carried out 
with different values of the cutoff radius, r c , for the 
nonbond interactions of the employed force field. The 
cutoff r c was changed from 11 A (the value for realistic 
conditions) to 4 A. Simulations were also performed 
for the isolated chain. The results reported in Ref. 
[87] are summarized in Fig. 29. Diminishing the cut- 
off radius yields a systematic decrease of the average 
barrier and the width of the barrier distribution g(Vs). 
An abrupt change of the latter parameters takes place 
in a narrow region around r c w 5 A, which roughly 
corresponds to the inter-macromolecular distances in 
PI. This result suggests that only the nonbond inter- 
actions at the inter-chain distances play a significant 
role in determining the width of the distribution of ro- 
tational barriers for methyl group rotation. Moreover, 
it is worthy of remark that these interactions also af- 
fect significantly to the average value of the potential 
barrier, i.e., the latter is not only determined by the 
chemical structure of the monomcric unit. 



VIII. CONCLUSIONS AND OUTLOOK 

This paper gives a comprehensive review of neutron 
scattering investigations on methyl group dynamics in 
glassy polymers of different chemical composition and 
microstructure. Considerable progress on this topic 
has been achieved over last 10 years mainly fueled by 
two facts: (i) the proposal in 1994 of the Rotation 
Rate Distribution Model for methyl group dynamics in 
polymers; (ii) the first direct experimental observation 
in 1998 by neutron scattering of quantum rotational 
tunnelling in glassy polymers at very low temperature. 

The experimental data presented in this paper 
demonstrate the uniqueness, the power, and the effi- 
ciency of the combination of quasielastic and inelastic 
neutron scattering for investigating this problem. In 
particular, it is worthy of remark that neutron scat- 
tering is the only technique which provides a direct 
experimental observation of quantum tunnelling pro- 
cesses of methyl groups. There are, however, some 
limitations, mainly related to the current resolution of 
the available high-resolution backscattering spectrom- 
eters (w 0.3 /xeV), which prevents the observation of 
clear rotational tunnelling features at very low temper- 
ature in polymers when the average rotational barrier 
is higher than w 1000 K. 

On the other hand, in this article it is shown that 
distribution effects are essential to understand methyl 
group dynamics — from quantum tunnelling to clas- 
sical hopping — in polymeric materials and in other 
highly disordered systems as well. The Rotation Rate 
Distribution Model (RRDM), which introduces a dis- 
tribution of potential barriers for methyl group rota- 
tion, provides a consistent description for the exper- 
imental features observed at different temperatures, 
including the crossover from quantum to classical be- 
havior and the librational energies as well. 

Taking into account the RRDM results on different 
polymer systems as well as on low molecular weight 
glasses, it is now evident that the distribution of po- 
tential barriers is determined by the intermolecular 
disorder and that the chain connectivity does not play 
any significant role. However, appart from these qual- 
itative results, it is clear that the actual origin of the 
distribution of potential barriers and its possible re- 
lationship with the structure factor characterizing the 
packing conditions, is one of the main points that still 
needs to be clarified. Concerning this question, fully 
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atomistic molecular dynamics simulations can play an 
important role in the near future. The results avail- 
able by now suggest that the width of the distribu- 
tion of potential barriers is mainly determined by the 
nonbond interactions at the inter-chain distances, at 
least in the case of polyisoprcnc. These results need to 
be confirmed by similar simulations in other realistic 
model systems as for instance, PMMA. 

Another way of addressing this question experimen- 
tally would be to carry out neutron scattering mea- 
surements on methyl group dynamics in polymers un- 
der confinement conditions, and to compare them with 
the corresponding results for the bulk state. To the 
best of our knowledge, a quantitative determination 
of the barrier distribution for methyl group rotation 
has never been reported for a confined polymer. How- 
ever, this has been recently carried out for glassy 
toluene confined in mcsoporous silicates [94]. From 
the knowledge of the barrier distribution in the bulk 
glassy state [35], and by investigating the dependence 
of the RRDM parameters on the monodisperse radius 
of the confining pores, a lower limit has been esti- 
mated for the distance beyond which methyl groups 
do not feel the interaction with the pore wall and 
behave as bulk-like rotors. This procedure has con- 
firmed that the interaction range relevant for methyl 
group dymamics in glassy toluene extends beyond the 
typical nearest-neighbour distance between centers-of- 
mass [94] . Similar investigations on confined polymers 
would be helpful to improve our understanding of the 
structural origin of the barrier distribution for methyl 
group rotation in polymers. 
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FIG. 1. Picture of methyl group rotation about its 
C3-axis, and level scheme for a purely threefold potential 
with barrier height Vz = 500 K. 



Detector 




FIG. 2. Scheme of a scattering experiment. The sphere 
of radius fco represents elastic interactions. Indexes 1 and 2 
respectively correspond to inelastic processes with neutron 
energy loss and gain. 




FIG. 3. Thick lines: A(Q) (solid) and l~A{Q) (dashed) 
for methyl group classical hopping. Thin lines: A(Q) 
(solid) and (1 — A(Q))/2 (dashed) for methyl group ro- 
tational tunnelling. 
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FIG. 4. Experimental spectra (circles) for methyl group 
dynamics in crystalline sodium acetate trihydrate. Solid 
lines are fits to the model for crystalline systems exposed 
in Section III. Dashed lines correspond to the experimental 
resolution. Q=1.8A . Scales (referred to the maximum): 
(a) 8%, (b)-(d) 5%. 



FIG. 5. Temperature dependence of the shift of the 
tunnelling frequency and of the crossover and classical 
Lorentzian HWHM for crystalline sodium acetate trihy- 
drate. Solid lines in panels (a-d) are fits to respectively, 
Eqs.(18),(23),(24) and (15). Panel (e) shows all the data 
in a common frame. The arrow indicates the crossover 
temperature T c . Dotted lines are a guide for the eyes. 
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FIG. 7. Solid lines: dependence on the barrier height 
of the tunnelling frequency (a) , first librational energy (b) 
and classical activation energy (c) for a protonated methyl 
group in a purely threefold potential. The dashed line in 
panel (a) corresponds to a deuterated methyl group. 
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FIG. 8. Distributions g(V 3 ), f{E A ), F(E 01 ) and h{hw t ) 
for an average barrier (V3) = 800 K and a standard devi- 
ation o y = 250 K. 
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FIG. 12. As Fig. 11, but KWW fits are done in different 
narrow intervals covering, as in usual experimental condi- 
tions, 1-2 decades in time. The obtained values of and 
to for the KWW description are shown. 
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FIG. 13. IN16 spectra for PVAc at Q = 1.8k 1 and dif- 
ferent temperatures. Circles are experimental data. Solid 
lines correspond to the RRDM description. Dashed lines 
are the resolution function. The scale is a 5% of the max- 
imum. 
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FIG. 15. Temperature dependence of To (a) and a (b) 
for PVAc. Different symbols indicate the different used 
spectrometers: IN16 (squares), IN5 (circles) and IN6 (tri- 
angles). Solid lines for To and a are respectively fits to Eqs. 
(45) and (46). The inset in panel (b) shows the values of 
<je obtained at each temperature from the corresponding 
values of o(T) by transformation through Eq. (46). The 
horizontal line corresponds to the value obtained by fitting 
the latter values to Eq. (46). Dashed lines in panel (a) 
correspond to Y = Too exp[— {(Ea) ± ffE)/feT]. 




FIG. 16. Generalized VDOS (circles) for PVAc, mea- 
sured at TOSCA at T = 20 K. The curve corresponds 
to the RRDM theoretical distribution F(Eoi). The arrow 
indicates the maximum of the latter for comparison with 
experiment. 




FIG. 17. Experimental integrated inelastic intensities 
for PVAc at Q = 1.8 A -1 in the intervals 1-2.5 fieV 
(squares), 2.5-6 /ieV (triangles) and 6-10 fieV (circles). 
Solid lines correspond to the theoretical description by the 
general version of the RRDM. Dashed lines are an extrap- 
olation of the classical limit of the RRDM to low temper- 
atures. 
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FIG. 18. Temperature dependence, for several poly- 
mers, of \n[I FE w(Q,T)/I FE w(Q,T « 0)], with 
Ifew(Q, T) the FEW intensity, measured in BS spectrom- 
eters. For clarity, data are consecutively shifted by 0.1. 
Q = 1.8 A' 1 for PDMS and PSF; Q = 1.95 A" 1 for PI 
and atactic PP. 
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FIG. 20. Methyl group incoherent contribution, 
(Q,lu), normalized to unity, for PMMA-d5 (circles) 
and PMMA-d8 (triangles) at T = 2 K and Q = 1.8 A -1 . 
Dashed lines correspond to the description by the tun- 
nelling limit of the RRDM. Solid lines correspond to the 
resolution function. Data for PMMA-d5 are shifted by 0.02 
for clarity. 




FIG. 19. Distribution of tunnelling frequencies for 
PMMA-d5 (dashed line), and PMMA-d8 (solid line). The 
inset shows both distributions in logarithmic scale. 
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FIG. 21. Distribution of tunnelling frequencies for sev- 
eral polymers, as derived from the RRDM analysis of hop- 
ping and/or librational data. 
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FIG. 22. Circles: Experimental spectrum at T — 2 K 
for PDMS. Thick and thin solid lines correspond to the 
theoretical spectra calculated respectively from the RRDM 
parameters reported in Refs. [89] and [93]. The dashed 
line is the resolution function. The scale is a 3 % of the 
maximum. 
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FIG. 23. Generalized VDOS (circles) measured at 
TOSCA, and RRDM theoretical distributions F(E i) 
(solid lines). All data are obtained at T = 20 K. (a): 
PDMS (b): PMMA-d5. Arrows indicate the maxima of 
the theoretical F(Eoi) for comparison with experiment. 
In panel (b), the thin line corresponds to a purely three- 
fold approximation, while the thick line corresponds to a 
sixfold correction Vq/Vs = 0.11 (see text). 
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FIG. 24. Generalized VDOS for PI measured at IN6. 
T = 140 K. Circles: PI-d5 (deuterated main chain - proto- 
nated methyl group). Triangles: PI-d3 (protonated main 
chain - deuterated methyl group). Curve: RRDM theoreti- 
cal distribution F(Eoi)- The arrow indicates the maximum 
of the latter for comparison with experiment. 
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FIG. 26. Symbols: Experimental spectra for a same 
sample of sodium acetate trihydrate in the crystalline (tri- 
angles) and glassy (circles) state at T — 2 K and Q = 1.8 
A -1 . The thick solid line is the RRDM theoretical func- 
tion. The thin solid line is a guide for the eyes. The dashed 
line is the resolution function. Scale: 3% of the maximum. 
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FIG. 25. Generalized VDOS for PSF measured at 
TOSCA. T = 20 K. Circles: fully protonated sample. Tri- 
angles: deuterated methyl group. Curves: RRDM theo- 
retical distribution F(Eoi) for a single-particle approach 
(dashed) and for an additional methyl-methyl coupling in- 
teraction (solid, see text). 
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FIG. 27. Circles: Experimental tunnelling spectrum at 
T = 2 K and Q = 1.8 A" 1 for tri(vinyl acetate). Solid line: 
Theoretical spectrum calculated from the RRDM parame- 
ters for PVAc. The dashed line is the resolution function. 
Scale: 3% of the maximum. 
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FIG. 29. Panel (a): Effect of the cutoff radius, r c , for 
the nonbond interactions on the distribution of potential 
barriers for methyl group reorientation in PI. Panel (b) 
shows the evolution with r c of the average potential barrier 
(V3} and the standard deviation cry. The shadowed area 
corresponds to the range of inter-chain distances. 
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TABLE I. Scattering cross-sections of the most abun- 
dant nuclei in polymers. Values are given in barn (1 barn 
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(E A ) (K) o e (K) (meV) (V z ) (K) a v (K) EST (meV) 



PVME [25,115] 


1010 


276 


7.8 










PVAc [31,33,89,146] 


445 


253 


9.1 


534 




274 


14.4 


PI [92,147] 


1170 


353 


62.1 


1364 




376 


23.4 


PDMS 

[89,146] 
[93] (*) 


842 
602 


144 
132 


18.5 
0.63 


997 




155 


21.2 


PEP [111] 


1288 


162 


1.6 










PP 

atactic [90,145] 
isotactic [144] 
head-to-head [111,145] 


1744 
1242 (§) 


421 
197 (§) 


28.8 
1.4 


1995 




448 


30.0 
28.5 
28.0 


DMB [124] 


842 


361 


40 










DHMS [124] 


1071 


385 


40 










PSF [149] 


2553 


396 


4.1 


2727 ; W£ = 


522 


406 


36.7, 41.6 (f) 


PH [149] 


2553 


373 


8.3 


2727 ; W£ = 


522 


383 


36.7, 41.6 (f) 


PC [149] 


2553 


379 


10.4 


2727 ; W£ = 


522 


390 


36.7, 41.6 (f) 


PMMA (ester-methyl) 
syndiotactic [32,34] 
SI: (s:i:h) = (50:10:40) [91] 
blend isotactic/Sl 
non-sc (i-methyl) [91] 
sc (i-methyl) [91] 
sc (Sl-methyl) [91] 


710 

529 

770 

686 
553 


241 
192 

313 
205 
180 


4.8 
1.5 

4.8 
2.8 
1.5 


820 ; (V 6 ) = 


90 


265 


22.4 


PMMA (a-methyl) 
isotactic [17] 
syndiotactic [17] 














37.2 
44.6 


PaMS 
syndiotactic [17] 

atactic [17] 
head-to-head [17] 














47.1 
47.1 
37.2 


PPO [16,128] 


1684 (t) 












28.3 


PIB [145,147,148] 






44 








40.0 



blend PVME/PS 
(65:35) [115] 
(20:80) [115] 



1010 
1010 



278 
354 



7.8 
7.8 



blend hh-PP/PEP 
(50:50; hh-PP-methyl) [111] 



1242 



197 



1.4 



n- hexat r iacont ane 
Mon (single-layered) [150] 
Orth II (double-layered) [150] 



30.0 
29.4 



TABLE II. Values of the RRDM parameters for the polymers investigated so far by neutron scattering techniques. 
Ea, cfe and Too are obtained from the analysis of high temperature data. Typical error bars are about a 10 % for the 
two formers and a 30 % for the latter. Eol ax denotes the maximum of the experimental librational peak. The parameters 
(V3) and ay of the barrier distribution g{Vz) are given in the cases where the quantities Ea, oe and Eo{ ax are known, 
and a good correspondence between them is possible by assuming a purely threefold potential. In the cases of PVAc and 
syndiotactic (s) PMMA also the tunnelling and crossover regimes have been investigated. For PVAc a good description 
of the latters is possible with the given RRDM parameters, but for s-PMMA sixfold corrections with ratio Ve/Vs = 0.11 
are necessary (see Section VI). (Ve) is indicated. In the case of PSF, PH and PC, a non-distributed coupling term is 
necessary to reproduce the librational double-peak (see Section VI). The coupling term W c is indicated. For the systems 
where hopping or librational data are absent, there is no sufficient information to determine the shape of the potential, 
and the parameters of g{Vs) are not given. 

(*): These parameteres, obtained from an analysis of hopping data, are uncompatible with both the experimental 
librational peak and the tunnelling spectrum (see Section VI). 

(§): In a purely threefold approach, these parameters provide a theoretical librational maximum of ~ 26 meV, close but 
not fully satisfactory with the experimental value (see Section VII). 
(f): Double-peak structure. 

(X)'- The standard deviation oe is not reported. In a threefold approach, the corresponding librational energy for the 
given (Ea) is 29.7 meV, to be compared with the experimental maximum. 
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